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ABSTRACT 

In  this  thesis  a  profitability  model  is  developed  for 
an  exploratory  well  drilling  program.   The  model  considers  the 
number  of  wells  drilled,  the  probability  of  success,  the  reserve 
size  discovered  by  a  successful  well,  the  discounted  profit  per 
barrel  for  oil  produced  and  the  cost  of  an  unsuccessful  well  or 
"dry  hole".   Each  of  these  parameters  are  analyzed  and  a  prob- 
ability distribution  developed  for  all  but  the  dry  hole  costs. 
With  this  model  and  the  respective  probability  distributions, 
a  FORTRAN  IV  computer  program  is  developed  using  the  Monte 
Carlo  technique  of  stochastic  simulation.   The  program  provides 
the  best  estimate  of  discounted  profit  for  the  venture,  the 
risk  of  the  venture  as  measured  by  the  standard  deviation  (and 
the  coefficient  of  variation)  and  the  probability  distribution 
of  the  discounted  profit  -  all  after  any  number  of  prespecified 
model  executions. 

After  developing  a  hypothetical  set  of  base  parameter 
data  and  executing  the  program,  a  reference  set  of  values  for 
the  profitability  and  risk  are  obtained.   Each  parameter  set 
is  then  varied  holding  the  remaining  parameters  at  their  base 
values.   Using  this  technique  the  effects  of  parameter  changes 
on  observed  model  results  are  assessed. 

Several  different  sensitivity  studies  are  performed  in 
this  manner  using  the  data  thus  developed.   These  studies 
include  the  effects  of  parameter  changes  on  the  expected  profit 
and  risk  of  the  venture  using,  a  Significance  Ratio  technique 
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developed  in  the  thesis.   In  addition,  the  convergence  of  the 
best  estimate  of  discounted  profit  with  the  number  of  model 
executions,  as  well  as  the  relationship  between  the  number  of 
wells  drilled  and  the  venture  risk  are  analyzed. 

Also  studied  is  the  possible  effects  of  periodicities 
in  the  combined  distribution  data  with  an  increased  number  of 
model  executions  and  the  effects  of  changing  a  parameter  prob- 
ability distribution  on  the  best  estimate  of  discounted  profit 
for  the  venture. 
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CHAPTER  I 
INTRODUCTION 

1 . 1   General 

A  primary  consideration  of  management  in  making  invest- 
ment decisions  is  the  eventual  profitability  of  the  proposed 
venture.   There  are  several  ways  of  measuring  profitability 
such  as  Rate  of  Return,  and  Present  Worth  or  Discounted  Expect- 
ed Profit.   The  main  purpose  in  calculating  this  profitability 
criterion  is  to  summarize  into  a  single  measure  the  quantifi- 
able factors  affecting  the  economic  desirability  of  the  project 
under  question  [25] . 

Profitability  calculations  have  generally  been  based 
on  the  best  estimates  of  the  individual  variables  involved. 
However,  these  best  estimates  may  or  may  not  occur.   The  re- 
sultant single  value  of  profitability  based  on  these  best  es- 
timates gives  no  information  on  the  chances  of  doing  much 
better  or  much  worse  then  the  expected  value. 

In  his  article,  Risk  Analysis  in  Capital  Investment 
[23],  David  B.  Hertz  demonstrates  that  in  certain  cases  the 
conclusions  reached  on  investment  decisions  may  have  been  much 
different  if  the  analysis  provided  not  only  a  single  profit- 
ability estimate  but  also  the  chances  of  that  estimate  being 
exceeded  or,  on  the  other  hand,  not  being  met.   Hertz  concludes 
that  the  application  of  probabilities  will  often  yield  entirely 
different  and  better  decisions.   Thus,  by  developing  the  pro- 
bability distribution  of  the  profitability,  management  can 


make  quantitative  assessments  of  the  risks  involved  in  approving 
a  particular  investment.   The  chances  of  achieving  a  minimum 
desired  return  or  of  losing  money  can  be  determined  so  that 
management  knows  the  size  of  the  risk  it  is  undertaking. 

In  order  to  determine  the  profitability  of  the  venture 
and  the  probability  distribution  of  the  profitability,  a  model 
must  be  established  containing  all  the  factors  which  bear  on  the 
problem.   In  some  cases  the  model  is  of  such  a  nature  that  it 
can  be  solved  analytically.   That  is,  the  model  consists  of 
known  mathematical  procedures  or  algorithms,  which  when  manip- 
ulated, either  directly  or  iteratively,  provide  a  solution. 
Many  types  of  managerial  problems,  however,  are  so  complex  that 
neither  direct  nor  iterative  analytical  procedures   exist  for 
their  solution,  or,  if  such  methods  do  exist,  they  represent 
such  great  difficulties  in  application  that  their  utilization 
is  impractical.   In  such  cases  we  turn  to  the  use  of  a  technique 
called  simulation.   Richards  and  Greenlaw  [42]  discuss  simula- 
tion as  follows: 

Simulation  involves  designing  and  utilizing  a  model 
which  replicates  some  aspect  of  the  firm's  operations... 
In  general,  simulation  models  may  be  distinguished  from 
other  operations  research  models  in  that  they  are  non- 
analytical  in  nature,  calling  for  an  experimental  app- 
roach, in  which  workable,  satisficing  solutions  rather 
than  optimal  ones  are  usually  provided. 

In  this  thesis  an  oil  exploration  investment  decision  model  will 
be  developed  incorporating  several  different  variables,  all  of 
which  bear  on  the  profitability  of  an  oil  exploration  venture. 


Each  variable  will  have  a  different  probability  distribution 
making  the  model  so  complex  that  analytical  solution  is  not 
practical.   We  will  therefore  turn  to  simulation  techniques, 
namely  the  Monte  Carlo  technique  (described  in  Chapter  II) ,  to 
provide  a  measure  of  profitability  and  a  probability  distribu- 
tion of  that  profitability. 

Before  discussing  the  organization  of  this  thesis  it 
may  be  well  to  clarify  one  point  which  seems  to  cause  a  great 
deal  of  confusion,  that  is,  the  distinction  between  uncertainty 
and  risk.   Contemporary  literature,  by  and  large,  still  maintains 
F.H.  Knight's  distinction  between  these  two  terms  [29].   Risk 
is  characterized  by  situation  in  which  the  outcome  is  not  known 
with  certainty,  but  where  the  probabilities  of  the  outcomes  are 
either  known  or  can  be  estimated.   Uncertainty,  on  the  other 
hand,  refers  to  situations  for  which  probabilities  of  the  out- 
comes cannot  even  be  predicted  in  probabilistic  terms. 

Concerning  this  distinction,  Alfred  Rappaport  [39], 

writes : 

While  this  distinction  may  be  a  useful  pedagogical 
device,  there  is  serious  doubt  about  its  applicability 
to  decision  makers  in  an  organizational  context.   Here 
the  distinction  is  blurred  by  the  fact  that  decision 
makers  generally  have  some  feelings  about  the  probabil- 
ities of  future  events.   Admittedly,  these  feelings 
range  from  a  degree  of  confidence  bordering  on  certainty 
to  ill-defined  feelings  approaching  mysticism.   The 
critical  point  remains  that  these  feelings  affect  the 
judgments  made  by  the  decision  maker. 

Based  on  the  reasons  expressed  by  Rappaport  the  terms  "uncert- 
ainty" and  "risk"  will  be  used  interchangeably  throughout  this 
thesis . 


1 . 2  Organization  of  the  Thesis 

This  thesis  treats  the  sensitivity  studies  in  profit- 
ability analysis  in  two  parts.   The  first  part,  beginning  with 
Chapter  II,  describes  the  Monte  Carlo  method  of  stochastic 
simulation,  a  technique  used  throughout  the  thesis.   In  Chapter 
III  a  general  profitability  model  for  an  exploratory  oil 
drilling  venture  is  developed.   Chapter  IV  discribes  some 
general  statistical  "yardsticks'*  used  to  measure  central  ten- 
dency and  variability  of  probability  distributions.   Beginning 
with  Chapter  V  and  continuing  through  Chapter  VIII  each  varia- 
ble, or  parameter,  in  the  profitability  model  is  analyzed  and 
a  probability  distribution  best  describing  the  parameter  is 
established.   In  addition,  a  hypothetical  set  of  values  for 
each  parameter  is  established  which  will  be  used  as  a  Base  for 
the  subsequent  sensitivity  studies.   With  the  profitability 
model  and  the  respective  probability  distributions  established, 
a  FORTRAN  IV  computer  program  is  developed  in  Chapter  IX  using 
the  Monte  Carlo  method  to  provide  the  best  estimate  of  dis- 
counted profit  for  the  venture,  the  probability  distribution  of 
the  profit  and  the  risk  of  the  venture.   Finally,  in  Chapter  X, 
using  the  hypothetical  set  of  values  for  each  parameter  pre- 
viously developed,  the  model  is  executed  and  a  reference,  or 
Base,  set  of  values  are  obtained  for  the  venture  profitability 
and  risk. 

The  second  part  of  the  thesis  deals  with  the  actual 
sensitivity  studies  using  the  model,  parameter  probability 


distributions  and  Base  values  developed  in  previous  chapters. 

Chapter  XI  describes  the  need  for  sensitivity  analysis 
in  general  terms  and  sets  the  framework  for  the  studies  which 
follow.   In  Chapter  XII  the  results  of  the  sensitivity  studies 
are  reported  and  analyzed.   Conclusions  and  recommendations 
resulting  from  these  studies  and  analyses  are  contained  in 
Chapters  XIII  and  XIV,  respectively. 


CHAPTER  II 
STOCHASTIC  SIMULATION  -  THE  MONTE  CARLO  TECHNIQUE 

2 . 1  History 

Through  the  ages  man  has  used  chance  processes. 
Theories  of  chance  processes  were  known  in  ancient  times.   The 
evidence  is  in  the  casting  of  the  astragal  bones  from  the  ankles 
of  sheep  and  goats  in  the  manner  of  dice  to  subdivide  land 
holdings.   This  was  called  casting  lots.   However,  it  was  not 
until  the  17th  century  that  probability  theory  became  formalized. 
Interestingly,  this  mathematical  development  arose  out  of  in- 
quiries into  games  of  chance. 

Another  important  development  occurred  in  the  1940' s 
during  war  work  on  the  atomic  bomb.   Scientists  were  frustrated 
by  being  unable  to  solve  some  complex  nonprobabilistic  mathe- 
matical problems  directly  by  analysis.   By  combining  chance 
processes  and  probability  theory,  they  solved  the  problems  in- 
directly.  To  do  this,  they  created  a  new  problem  stated  in  terms 
of  stochastic  models  in  such  a  way  that  the  new  problem  was  equ- 
ivalent to  the  complex  nonprobabilistic  problem.   But  the  sto- 
chastic version  represented  chance  processes  that  could  actually 
be  carried  out.   They  actually  performed  the  chance  processes 
represented  by  the  probabilistic,  or  stochastic,  formulas  in 
the  new  version  of  their  problem.   Their  methods  were  equivalent 
to  flipping  coins  or  throwing  dice,  but  much  more  efficient. 
The  results  of  these  experiments  were  then  substituted  for  the 
stochastic  expressions  of  the  new  version  of  the  problem  to  give 
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answers  to  the  original  nonprobabilistic  problem.   These  scien- 
tists gave  the  name  Monte  Carlo  to  their  new  technique  [5]  [22] . 

2 . 2  The  Monte  Carlo  Technique 

Monte  Carlo  techniques  were  born  when  probability  theory 
and  actual  chance  processes  were  used  to  solve  problems  that  had 
no  stochastic  aspects  at  all.   Since  that  time  Monte  Carlo  tech- 
niques have  been  used  to  solve  a  multitude  of  nonstochastic  (or 
deterministic)  type  problems  [33] . 

When  dealing  with  stochastic  (or  probabilistic)  simula- 
tion it  is  necessary  to  have  a  means  for  dealing  with  those 
probabilistic  variables,  the  values  which  assume  a  frequency 
distribution.   Richard  F.  Barton  [4]  defines  stochastic  simu- 
lation as: 

. . .  one  in  which  differing  outputs  trial  to  trial  can 
be  obtained  without  changing  the  inputs  (ignoring  random 
numbers  as  inputs) .   Specifically,  this  means  that 
identical  parameters,  starting  conditions,  and  input 
time  path  values  produce  varying  outputs  trial  to  trial 
and  run  to  run. 

One  commonly  used  method  for  working  with  stochastic  variables 

in  simulation  models  is  again  the  Monte  Carlo  technique.   The 

technique  is  one  which  both  has  application  to  an  almost  endless 

variety  of  stochastic  simulation  problems  and  at  the  same  time, 

is  relatively  simple  to  comprehend  and  easy  to  utilize.   S.W. 

Hess  and  H.A.  Quigley  [24]  define  the  Monte  Carlo  technique  as: 

...  a  sampling  procedure  whereby  complicated  expressions 
involving  one  or  more  probability  distributions  may  be 
evaluated.   In  essence,  it  consists  of  simulation  an 
experiment  to  determine  some  probabilistic  property  of 
a  population  of  objects  or  events  by  the  use  of  random 
sampling  applied  to  the  components  of  the  objects  or 
events . 


A  parameter  value  is  obtained  by  drawing  randomly  from 
the  parameter's  probability  distribution.   In  a  similar 
manner,  selections  are  made  for  each  parameter  value 
from  its  respective  distribution.   This  set  of  parameter 
values  is  then  substituted  into  the  model  and  the  first 
sample  value  of  the  independent  variable  is  computed. 
Subsequent  values  of  the  independent  variable  are  obtain 
ed  by  drawing  additional  sets  of  parameter  values. 

The  individual  values  of  the  independent  variable  so 
obtained  are  in  approximation  to  its  true  probability 
distribution.   The  approximate  distribution  approaches 
the  true  distribution  as  the  sample  size  increases. 

The  Monte  Carlo  technique  actually  involves  a  two  sta- 
ge process.   In  the  first  stage,  stochastic  relationships  be- 
tween the  variables  are  determined  in  order  to  construct  the 
simulation  model.   These  relationships  may  be  determined  from 
surveys,  company  data,  industry  data,  executive  experience  and 
the  like.   At  any  rate,  once  the  relationships  are  determined, 
probabilities  of  occurrence  of  frequency  distributions  are 
assigned  for  each  variable.   The  frequency  distributions  are 
then  represented  by  associating  groups  of  numbers  to  each  value 
of  the  distribution,  the  size  of  the  group  being  governed  by 
the  relative  frequency  of  the  value.   As  an  example,   assume 
60%  of  the  items  in  the  distribution  have  a  value  A,  30%  have 
a  value  B  and  10%  have  a  value  C.   These  three  values  can  be 
represented  by  numbers  between  00  and  99.   Thus  the  items 
having  a  value  A  would  be  represented  by  the  numbers  from  00 
through  59,  those  having  values  of  B  and  C  would  be  represented 
by  numbers  from  60  through  89  and  90  through  99,  respectively. 

Once  the  model  is  developed,  frequency  distributions 
assigned  and  numbers  representing  these  frequency  distributions 


determined,  we  move  to  stage  two  which  involves  running  the 
Monte  Carlo  process  to  obtain  a  variable  value.   We  do  this  by 
first  selecting  a  random  number  (or  numbers)  from  tables  if  the 
simulation  is  to  be  performed  manually  or  by  building  a  random 
number  generator  into  a  computerized  simulation  program. 
(With  some  computer  systems,  such  as  the  CDC  6600,  it  is  poss- 
ible to  make  use  of  a  random  number  generating  function  in- 
corporated in  the  computer  systen  itself.)   Once  the  random 
number  is  selected  (or  generated)  it  is  matched  to  the  corres- 
ponding number  assigned  to  represent  each  value  of  the  frequency 
distribution  of  the  variable  in  question.   This  in  turn  estab- 
lishes the  value  of  that  variable.   This  process  is  repeated 
for  each  variable  in  the  model,  the  variables  combined  and  a 
single  answer  is  obtained.   Stage  two  is  then  repeated  as  many 
times  as  is  necessary  to  obtain  a  statistically  reliable  sample. 
The  various  stages  of  the  Monte  Carlo  process  have  been  schem- 
atically diagrammed  by  M.D.  Richards  and  P.S.Greenlaw  [41] 
and  are  shown  in  Figure  2.1. 
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CHAPTER  III 
THE  PROFITABILITY  MODEL 

3 . 1  General 

As  noted  in  Chapter  I,  the  profitability  of  an  explor- 
atory well  drilling  venture  and  the  probability  distribution 
of  that  profitability  will  be  analyzed  using  a  model  which 
incorporates  the  several  different  factors  which  bear  on  the 
problem,  namely:   the  number  of  wells  drilled  in  the  program, 
the  probability  of  a  success,  the  size  of  field  discovered  by 
a  successful  well,  the  discounted  profit  per  barrel  of  oil 
discovered  and  the  cost  of  a  dry  hole,  or  an  unsuccessful  well. 
The  model  essentially  measures  the  profitability  of  the  venture 
on  the  basis  of  Net  Present  Worth  technique  whereby  the  expect- 
ed cash  flows  are  discounted  at  a  suitable  discount  rate  to 
reflect  the  time  value  of  money.   In  the  model,  and  throughout 
this  thesis,  the  following  notation  will  be  utilized: 

P  =  Total  discounted  profit  on  any  one  drilling  endeavor. 

E(P)  =  The  best  estimate  (mean)  of  discounted  profit  for 
the  total  exploratory  well  drilling  venture. 

n  =  Number  of  wells  drilled  in  the  drilling  venture. 

p  =  Probability  of  a  successful  well. 

x  =  Number  of  successful  wells,  where  x  can  assume  values 
from  o  to  n. 

Q  =  Reserve  size  of  oil  field  discovered  by  a  successful 
well  (in  million  bbls) . 
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ZQ  =  Reserve  size  of  oil  fields  discovered  in  the  venture 

(in  million  bbls) . 
e  =  Discounted  profit  per  barrel  for  oil  produced 
(in  $/bbl). 
Cu  =  Cost  for  an  unsuccessful  well,  or  dry  hole. 

Using  this  notation,  the  model  for  a  single  drilling  endeavor 
might  be  expressed  as  follows: 

P  =  x  Q  e  -  (n-x)  Cy  (1) 

In  an  exploratory  well  drilling  venture  it  must  be 
realized  that  the  number  of  successful  wells  is  a  function  of 
both  the  number  of  wells  drilled,  n,  and  the  probability  of 
success,  p.   Therefore: 

x  =  f(n,  p)  (2) 

Further,  the  size  of  reserve  discovered  in  the  venture, £  Q, 
is  a  function  of  the  number  of  successful  wells. 

EQ  =  f(x,  Q)  (3) 

It  will  be  assumed  that  discounted  profit  per  barrel,  e, 
and  dry  hole  cost  per  unsuccessful  well,  both  as  absolute 
values,  are  independent  in  that  they  are  not  directly  derived 
from  the  number  of  wells  drilled,  etc.   Therefore: 

P  =  [sQ  -  f  (x,  Q)]  e  -  [  n  -  [x  =  f  (n,  p)j]  Cu      (4) 

In  Equation  4,  which  will  be  used  as  the  profitability 
model  in  this  thesis,  it  can  be  seen  that  the  discounted  profit 
for  the  well  drilling  venture  is  equal  to  the  size  of  reserves 
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discovered  in  the  venture  (which  is  a  function  of  the  number 
of  successful  wells  and  the  reserve  discovered  by  each  success- 
ful well)  multiplied  by  the  discounted  profit  per  barrel  for 
oil  which  will  be  produced  from  the  reserves,  less  the  number 
or  unsuccessful  wells  (a  function  of  the  number  of  wells 
drilled  and  the  probability  of  successes)  multiplied  by  the 
cost  of  the  unsuccessful  well,  or  dry  hole  cost. 

In  terms  of  best  estimates,  the  best  estimate  of 
discounted  profit  for  the  venture,  E(P),  can  be  written, 

E(P)  =  E(x)  [E(Q)  E(e)]  -  [n  -  E(x)]  Cu  (5) 

The  relationship  expressed  in  Equation  5,  has  been  used 
by  many  to  estimate  the  profitability  of  a  venture  directly  as 
follows : 

Based  on  all  available  information,  "best  guesses"  are 
made  of  the  probability  of  success,  which,  in  this  case,  is 
the  chance  that  at  least  one  oil  field  will  be  found,  the 
magnitude  of  the  reserves  which  would  be  expected  by  a  success- 
ful well,  the  profit  per  barrel  for  oil  produced  and  the  dry 
hole  cost.   It  should  be  emphasized  that  normally  only  one 
value  for  each  parameter  is  provided,  that  is  the  best  estimate 
(or  guess) .   These  best  estimates  are  then  combined  with  the 
number  of  successes  determined  by 

x  =  np  (6) 

That  is,  if  we  had  an  18  well  program  and  the  probability  of 
success  was  1  in  9  then  we  could  expect  2  successes.   If  we 
multiply  the  2  successes  by  the  best  guess  of  reserve  size  per 
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well  (say  5  million  barrels)  we  would  have  a  total  reserve  of 
10  million  barrels.   If  we  estimate  the  most  likely  discounted 
profit  per  barrel  to  be  $.80  then  we  would  expect  10  million 
bbls  x  $.80/bbl  or  $8  million.   From  this  we  would  have  to 
deduct  dry  hole  costs  which  would  be  (18-2)  or  16  wells  at, 
say,  $600,000  per  well,  or  $9.6  million.   The  best  estimate  of 
discounted  profit  for  the  venture  based  on  this  technique  is 
therefore  -  $1.6  million.   In  all  probability  the  venture  would 
not  be  started  on  the  basis  of  this  analysis. 

As  discussed  in  Chapter  I,  what  this  analysis  ignores 
is  the  range  of  values  possible  for  each  of  the  parameters. 
For  example  the  best  estimate  for  reserve  size  was  5  million 
barrels  per  success.   That  may  be  the  best  estimate,  but  what 
if  the  reserves  actually  exceed  this  best  estimate  by  251? 
The  best  estimate  of  discounted  profit  would  then  have  been 
$0.8  million  and  a  potentially  profitable  venture  may  have 
been  passed  up. 

In  this  thesis  the  basic  relationship  expressed  in 
Equation  4  will  be  utilized  for  the  profitability  analysis. 
In  subsequent  Chapters  each  of  the  parameters  will  be  explored 
to  determine  their  functional  relationships  and  proper  frequency 
distributions. 


CHAPTER  IV 
DESCRIPTIVE  MEASURES  OF  PROBABILITY  DISTRIBUTIONS 

4.1   Measures  of  Central  Tendency   (Mean,  Median  and  Mode) 

In  analyzing  each  parameter  in  the  profitability  model 
as  well  as  the  resultant  discounted  profit,  it  is  necessary  to 
know  not  only  its  probability  distribution  but  also  some  single 
measure  typical  of  the  parameter.   This  single  measure  will  be 
used  to  describe  the  central  tendency,  or  the  general  region  of 
continuum  in  which  the  distribution  lies. 

The  most  common  measure  of  central  tendency  is  called 
the  Mean  or  Expected  Value  and  is  designated  x  (or  x  )  or  E(x). 
The  expected  value  is  a  measure  of  location  in  the  sense  that 
it  roughly  locates  a  middle  or  "average"  value  of  the  variable, 
x  [14] .   The  mean  is  the  sum  of  a  set  of  values  divided  by  the 
number  of  values  summed.   Expressed  mathematically  the  mean  of 
a  set  of  n  numbers,  x-^ ,  X2 ,  x-^,  ....,  xn  is: 

.    *1  +  *2  +  x3  +  x4  •••  +  xn        (1) 

E(x)  =  x  =  v  * 

K    J  n 

Or,  where  f^  is  the  frequency  of  each  observation  of  x  (i.e., 
x^ ,  X2,  x^,  etc.)  designated  as  x^  we  can  write: 

E(X)  =  x  =  ± (2) 

v  J  n 

A  second  measure  of  central  tendency  is  the  Median. 
The  median  of  a  set  of  values  is  the  middle  value  when  the 
observations  are  arranged  in  either  ascending  or  descending 
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order.   The  median  is  sometimes  defined  as  the  value  above  and 
below  which  an  equal  number  of  values  occur  [31].   On  a  cumu- 
lative probability  (or  frequency)  plot,  the  median  corresponds 
to  the  value  opposite  501  cumulative  probability  (or  frequency) . 
The  third  measure  of  central  tendency  is  the  mode .   The 
mode  is  defined  as  the  most  frequently  appearing  value.   It  is 
also  referred  to  as  the  "most  likely  value".   It  should  be  noted 
that  a  distribution  does  not  have  to  have  a  mode  (i.e.  the  uni- 
form distribution).   Conversely,  a  distribution  may  have  more 
than  one  mode  (hence  the  designated  "bimodal"  for  a  distribution 
with  two  modes) . 


Why  is  there  a  difference  between  these  three  measures 
of  central  tendency?   Lawrence  Mann[32]  describes  this  diff- 
erence as  follows: 

This  difference  is  due  to  the  fact  that  the  distribution 
is  somewhat  skewed.   Skewness  is  the  degree  to  which  a 
distribution  departs  from  symmetry.   If  a  distribution 
leans  to  the  right  then  it  is  said  to  be  skewed  to  the 
left  (negative  skewness);  if  it  leans  to  the  left,  then 
it  is  said  to  be  skewed  to  the  right  (positive  skewness) 
See  Figure  4-1. 
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In  this  thesis  we  will  analyze  the  parameters  of  the 
profitability  model  in  terms  of  their  mean,  median  and  mode, 
as  appropriate.   In  addition,  the  expected  value  and  median 
for  each  case,  or  parameter  set,  will  be  measured.   The  expected 
value  (mean)  will  be  considered  the  best  estimate  of  discounted 
profit  for  each  parameter  set  and  will  be  used,  along  with  the 
median,  to  compare  the  various  cases  in  the  sensitivity  studies 
which  follow. 

4.2   Measures  of  Variability   (Variance  and  Standard  Deviation) 
The  expected  value  of  a  random  variable  is  an  "average" 
value.   However,  this  measure  provides  no  information  about  the 
variability  of  the  values  of  the  variable.   A  measure  of  this 
variability,  or  the  "spread"  or  "dispersion"  of  the  values  of 
the  random  variable,  adds  another  dimension  to  the  description 
of  the  random  variable. 

The  most  useful  quantification  of  uncertainty  associ- 
ated with  a  future  event  is  the  standard  deviation,  a.   It 
measures  the  extent  of  the  spread  between  individual  possible 
outcomes  of  the  event  and  the  expected  value  of  the  event. 
Another  measure  of  variability  is  provided  by  the  variance , 
found  by  squaring  the  standard  deviation. 

For  a  discrete  random  variable,  x,  the  variance  is: 


a. 


2  =  E  (x-x)2   =  i 


n       .  7 


n-l  (3) 
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where  x  =  the  mean  value  of  x 
x-  =  any  value  of  x 

n  =  number  of  values  of  x 
f^  =  frequency  of  x^ 
The  denominator  in  Equation  2  is  replaced  by  n-1  when  the 
statistics  are  carried  out  on  a  sample  size  n  taken  from  a 
large  population  -  when  trying  to  draw  conclusions  concerning 

the  large  population.   Division  by  n-1  in  Equation  3  makes 

2 
the  sample  variance,  sx  ,  an  "unbiased"  estimate  of  the  pop- 

ulation  variance  aY  . 

If  x  is  a  continuous  variable,  the  variance  is: 

ax2  ■  E   (x-x)2  =   /   (x^x)2  g(x)  dx        (4) 

-  00 

Evaluating  the  usefulness  of  variance  as  a  measure  of  dispersion, 

Samuel   Goldberg  [15]  remarks: 

One  difficulty  with  the  variance  is  that  it  does  not 
measure  dispersion  in  the  same  units  as  the  values  of 
x.   Thus,  if  x  has  dollar  values,  then  E(x)  is  a  certain 
number  of  dollars,  but  since  the  variance  is  the  mean 
square  deviation,  Var  (x)  is  measured  in  dollars 
squared.   It  is  in  order  to  have  a  measure  of  disper- 
sion in  the  same  units  as  the  values  of  x  that  we 
define  the  standard  deviation  as  the  square  root  of 
the  variance. 

Hence,  variance  can  be  adjusted  for  this  disadvantage  by  taking 

its  square  root,  thus  defining  the  standard  deviation,  a.   This 

adjustment  makes  the  standard  deviation  directly  comparable 

with  values  of  the  random  variable.   The  standard  deviation 

takes  into  account,  as  does  the  expected  value,  not  only  the 

extreme  values  but  all  intermediate  values,  logically  relating 
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uncertainty  to  the  expected  value  [45] . 

In  this  thesis  the  standard  deviation  for  the  Base  and 
all  other  parameter  sets  is  computed  by  first  calculating  the 
variance  using  a  relationship  derived  from  Equation  3  and  then 
taking  its  square  root  as  follows: 


(5) 


4. 3  Coefficient  of  Variation 

A  knowledge  of  the  expected  value  and  standard  deviation 
gives  much  information  about  a  random  variable,  even  though  one 
does  not  have  any  knowledge  of  the  individual  values  of  a 
random  variable.   However,  standard  deviation  is  a  measure  of 
absolute  uncertainty,  that  is,  it  measures  the  actual  amount 
of  variation  present  in  the  random  variable;  and  it  is  depen- 
dent on  the  scale  of  measurement.   There  is  a  need  for  additional 
important  information  to  facilitate  comparison  of  variation  of 
several  random  variables  [44] . 

The  measure  which  relates  uncertainty  to  the  expected 
value  of  the  random  variable  is  the  coefficient  of  variation, 
cv,  the  formula  for  which  is  as  follows: 

ax 
cv  =  -£  (6) 

X 

where  cv  denotes  the  coefficient  of  variation.   Use  of  this 
measure  of  relative  uncertainty  makes  the  variability  figures 
directly  comparable.   It  describes  the  degree  of  uncertainty 
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or  relative  risk  of  the  individual  expectations  about  the  ex- 
pected value  of  the  random  variable  [38]  . 

Weston  and  Brigham  [47]  describe  the  use  of  the  coeff- 
icient of  variation  in  the  following  passage: 

. . .  Use  the  coefficient  of  variation  to  compare  the 
riskness  of  alternative  investments  whenever  the  stand- 
ard deviation  would  be  misleading.   In  general,  the 
standard  deviation  is  appropriate  when  percentage  rates 
of  return  are  considered,  while  the  coefficient  of 
variation  is  appropriate  if  dollar  returns  are  consid- 
ered.  The  reason  for  this  distinction  is  that  rates  of 
return  are  already  "standardized"  for  size  of  invest- 
ment, but  some  form  of  standardization  is  necessary 
when  dollar  returns  are  concerned. 

In  the  sensitivity  studies  which  follow,  the  risk 
involved  with  each  Case,  or  parameter  set,  will  be  measured 
in  terms  of  its  absolute  amount  (standard  deviation)  and  in 
terms  of  its  relative  value  (coefficient  of  variation) .   These 
values  will  then  be  related  to  the  standard  deviation  and 
coefficient  of  variation  of  the  Base  set  of  parameters  indicat- 
ing the  change  in  risk  associated  with  each  parameter  change. 


CHAPTER  V 
PROBABILITY  OF  SUCCESS 

5  . 1  General 

In  assigning  the  probability  of  occurrances  of  future 
events  one  normally  starts  by  examining  past  experiences  with 
similar  events. 

The  petroleum  industry  has  been  very  methodical  in 
recording  drilling  statistics  for  years.   These  statistics  will 
provide  the  decision  maker  with  the  necessary  data  to  compute 
probability  of  success,  or,  as  they  are  called  in  the  industry 
the  national  success  ratios.   These  national  success  ratios  are 
helpful  but  may  be  misleading.   They  include,  for  example,  all 
wells  that  produced  any  amount  of  oil  or  gas,  regardless  of 
whether  the  amount  was  small  or  large,  wells  drilled  in  different 
formations,  shallow  and  deep  wells,  and  wells  drilled  with 
various  kinds  of  skills. 

While  the  national  success  ratios  provide  one  set  of 
past  experiences  they  may  be  of  little  help  in  assuming  prob- 
abilities for  a  particular  venture  in  a  specific  location. 
For  that  particular  venture  local  experiences  rather  than 
national  experience  might  be  more  helpful.   But  even  with  local 
experience  there  are  still  problems  in  making  the  statistics 
selective  so  as  to  have  relevance  to  the  venture  at  hand.   For 
instance,  should  all  wells  in  all  geological  provinces  be  in- 
cluded?  Or  only  those  that  are  producing  from  the  target 
horizon  being  sought  in  the  venture  under  consideration? 
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These  are  difficult  questions  to  answer  and  may  be  complicated 
by  the  fact  that  if  we  wanted  to  exclude  all  wells  not  producing 
from  a  specific  geological  province  a  breakdown  of  data  per- 
mitting this  exclusion  might  not  be  available. 

It  should  be  noted  that  while  local  and  national  ratios 
of  successes  do  provide  guides,  the  assignment  of  probability 
of  success  for  a  particular  venture  is  largely  subjective,  and 
in  some  cases,  intuitive.   But,  no  matter  how  the  decision  is 
arrived  at,  a  probability  of  success  for  the  venture  must  be 
stated  if  a  profitability   analysis  is  to  be  made.   Before 
examining  how  the  probability  of  success  will  be  used  in  the 
profitability  analysis,  let  us  examine  some  concepts  associated 
with  probabilities  as  they  relate  to  business  decisions. 

For  a  long  time  only  interpretation  of  probabilities 
has  been  the  classical  long-run,  relative  frequency  arguement, 
i.e.,  repetition  of  the  event  over  and  over  under  identical 
conditions.   For  the  majority  of  business  problems,  therefore, 
this  concept  has  been  almost  useless  as  a  practical  guide  for 
action.   Most  business  decisions  concern  singular,  or  only 
occasional,  events  that  may  never  be  repeated  -  particularly 
under  identical  conditions.   Therefore,  choosing  an  act  because 
it  offers  the  best  "long-run  average"  fits  only  a  limited  number 
of  business  decisions.   Yet  the  "long-run"  philosophy  persists. 
During  his  research,  Grayson  [18]  reported  the  following 
comments  by  operators: 
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If  you  drill  enough  holes,  the  law  of  averages  will 
take  care  of  you.  If  you  keep  playing  the  averages 
the  percentages  will  hit  sooner  of  later. 

When  he  questioned  them  as  to  what  they  meant  by  the  words 

"sooner  or  later",  "average",  etc.  most  operators  could  not 

be  specific. 

If  we  say  that  the  probability  of  success  is  1  in  7 

does  this  mean  that  if  we  drill  7  wells  we  will  have  one 

success?   The  answer  is  no!  It  means  that  if  an  event, such  as 

drilling  a  well,  is  repeated  over  and  over,  the  frequency  of 

success  will  tend  to  approach  the  1  in  7  success  ratio.   It  is 

therefore  important  to  remember  that  the  number  of  successes 

is  dependent  on  the  number  of  wells  drilled  as  well  as  the 

probability  of  success.   If  we  consider  the  1  in  7  ratio  an 

average  success  ratio  then,  as  Grayson  writes  [20] 

...  an  operator  who  can  drill  only  a  small  number  of 
wells  may  expect  large  variations  from  the  "average". 
The  implications  for  a  man  with  limited  funds  are  clear 
If  he  had  unlimited  funds,  or  even  a  very  large  amount 
of  funds,  he  might  be  able  to  survive  these  variations 
early  losses  or  long  "runs"  of  failures  -  while  waiting 
for  the  "average  to  come  in".   But  with  limited  funds, 
an  operator  may  be  ruined  before  the  "average"  is 
realized.   It  would  be  similar  in  example  to  the  man 
who  drowned  while  crossing  a  stream  with  an  "average" 
depth  of  two  feet. 
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5 . 2  The  Binomial  Distribution 

An  event  that  has  a  certain  probability  of  occurrence, 
p,  in  a  definite  number  of  trials,  n,  is  said  to  be  binomially 
distributed. 

When  we  drill  an  exploratory  well  for  oil,  one  of  two 
outcomes  is  possible.   We  can  have  a  successful  well,  which  as 
we  will  see  in  the  next  section  is  defined  as  a  well  which 
discovers  a  field  of  1  million  barrels  or  greater,  or  we  can 
have  an  unsuccessful  well,  or  a'dry  hole'.   If  we  let  p  equal 
the  probability  of  success  than  1-p,  or  q  ,  is  the  probability 
of  a  failure.   When  two  exploratory  wells  are  drilled  three 
things  can  happen:   both  can  be  successful,  one  can  be  success 

ful  and  one  dry,  or  both  can  be  dry  holes.   The  chance  that 

9 

both  are  successful  is  pxp  or  p  ,  that  one  is  successful  and 

1  dry  ispxq+pxqor  2pq  and  that  both  are  dry  is  q  x  q 
or  q  .   Using  this  same  logic,  if  three  exploratory  wells  are 
drilled  we  can  have  3  successes  (p  ) ,  2  successes  and  one  dry 
hole  (3p2q) ,   1  success  and  2  dry  holes  (3pq2)  or  3  dry  holes 
(q3) .   Writing  these  terms  in  equation  format  we  have: 

p3  +  3p2q  +  3pq2  +  q3  =  (p  +  q)3  (1) 

Thus,  if  n  trials  are  made  each  of  which  results  in  either 
"success"  or  "failure",  if  the  probability  of  success  on  each 
trial  is  p  and  of  failure  on  each  trial  is  q=l-p,  and  if  the 
trials  are  independent,  then  the  probabilities  of  all  possible 
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results  of  such  a  trial  are  given  by  the  various  terms  of  the 
binomial  expansion. 

(p+q)n  =  pn  ♦  np""1,  *  a%llip"-2q2  ..     nl     x  n-x  ..  +qn>  (2) 

f  x i ^n-x  j ! 

Where  the  fourth  term  in  Equation  2  is  the  general  term  and 
describes  the  probability  of  x  successful  exploratory  wells 
in  n  independent  trials.   Hence,  the  binomial  distribution: 

n!    „x„n-x 


«W  "  xTTW%T  PA1""A  (3) 


n! 


The  term  x;  fn-x")  !  ^s  °ften  written  as  nCx  (the  number  of  com- 
binations) .    Equation  3  can  therefore  be  expressed  as: 


w(x)  =  nCx  pxqn"X  (4) 

Again  it  is  emphasized  that  Equation  3  provides  the 
probability  of  x  successes  in  n  trials  only  if  the  probability 
of  success,  p,  is  constant  from  trial  to  trial  and  the  trials 
are  independent. 

To  illustrate  the  use  of  the  binomial  distribution  in 
predicting  successful  events  let  us  examine  the  following 
example . 

Suppose  we  are  planning  an  exploratory  well  drilling 
program  of  20  wells  and  we  estimate  the  probability  of  success, 
p,  to  be  1  in  7  or  say  15%.   Then  we  can  also  say  that  the 
probability  of  failure,  q  is  1-p,  or  85%.   Using  Equation  4, 
we  can  write  the  probability  of  having  x  success  in  20  trials 
as : 
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20-x 


w(x)  =  20Cx(.15)x(.85)zu"x  (5) 


To  determine  the  probability  of  0  successes  we  would  set 
x=0,  and, 

w(0)  =  20C0(.15)°(.85)20"° 

w(0)  =  — 2JL! (.15)°(.85)20 

0!  (20)!       *"    J 

w(0)  =  (.85)20 

w(0)  =  .0388   or   3.88% 

Likewise,  the  probability  of  1  success  (x=l)  can  be  written 

19 


W(1)  =  1!  (20-1)!  CIS^C.SS) 

2.4329  x  1018  iq 

w(l)  =  TT  (^15)(.85)iy 

1.2165  x  10i/ 


w(l)  =  20(.15) (.0456) 

w(l)  =  .1368   or   13.68 

A  similar  analysis  could  be  performed  for  each  success  until  a 
total  of  20  were  calculated.  The  terms  could  however,  be  read 
directly  from  tables  [43] .  For  this  example  we  would  have  the 
following  probability  for  each  possible  outcome: 
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TABLE  5-1 
PROBABILITY  OF  SUCCESS  -  BINOMIAL  DISTRIBUTION 


Number  of 
Successes  (x)        w(x)  in  percent 


0  3.88 

1  13.68 

2  22.93 

3  24.28 

4  18.21 

5  10.28 

6  4.54 

7  1.60 

8  0.46 

9  0.11 

10  0.03 

11  0.00 


20  0.00 

5 . 3  Simulation  of  the  Binomial  Distribution 

In  the  Monte  Carlo  studies  to  follow  we  will  want  to 
simulate  the  binomial  distribution  in  order  to  determine  the 
number  of  successful  wells  for  each  drilling  venture.   For 
values  of  p  near  50%  and  for  large  n  it  can  be  shown  that  the 
Binomial  distribution  approaches  the  Normal  distribution  [12] 
The  mean  value,  xQ ,  and  standard  deviation,  a,  are  given  by: 

xQ  =  np  (6) 


o   =  /npq  (7) 
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If  xQ  corresponding  to  a  cumulative  frequency,  W(x)  =  50%, 
which  is  a  shorthand  writing  for  W(x  £  x . )  =  50%,  and  x  +  a 
corresponding  to  a  cumulative  frequency,  W(x)  =  161  are  plotted 
on  normal  probability  paper  and  a  straight  line  drawn  between 
them,  the  probabilities  associated  with  all  other  values  of  x 
could  be  read  directly  from  the  graph.   Alternatively,  the 
Normal  distribution  function 


t    ,       1  (x-x0) 

w(x)  =  exp  ' 


ax  /2?     (    2ax2  J  (8) 

could  be  solved  analytically  to  obtain  the  value  of  w(x)  for 
each  success,  x. 

If  the  value  of  p  varies  substantially  from  50%  and  if 
the  value  of  n  is  relatively  small  can  the  normal  distribution 
be  used  to  simulate  the  binomial  distribution?   To  answer  this 
question  we  will  continue  the  analysis  of  the  20  well  drilling 
venture  that  had  a  probability  of  success  of  15%.   From  Table 
5-1  a  cumulative  probability,  W(x) ,  can  be  calculated  such  that 
W(x  <   x^)  is  the  chance  that  the  number  of  successes  will  be 
x-  or  smaller. 
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TABLE  5-2 
CUMULATIVE  PROBABILITY  OF  SUCCESS  -  BINOMIAL  DISTRIBUTION 

Number  of  Successes  (x.)   w(x-)  in  Percent    (W(x  <  x-)in  Percent 

0  3.88  3.88 

1  13.68  17.56 

2  22.93  40.49 

3  24.28  64.77 

4  18.21  92.98 

5  10.28  93.26 

6  4.54  97.80 

7  1.60  99.40 

8  0.46  99.86 

9  0.11  99.97 

10  0.03  100.00 

11  0.00  100.00 


.   20  0.00  100.00 

The  values  of  (Wx  £  x.)  corresponding  to  values  of  x 

on  Table  5-2  were  plotted  on  normal  probability  paper  (Figure 

5-1)  as  circular  points.   If  these  values,  in  fact  approached 

a  normal  distribution  then  these  points  would  fall  along  a 

straight  line  passing  through  x  at  (Wx  <,  x.)  =  50%  and  x  +  o 

at  W(x  1  x.)  =  841. 
l 

From  Equation  (6) 
x0  =  np 
x0  =  20(.15) 

*o  =  3 
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From  Equation  (7) 

a  =  /  npq 

a  =  •  (20) ( . 15)  ( . 85) 


a  =  /  2.55 
a  =  1.60 
Therefore,  xQ  +  a  =  4.60 

A  straight  line  is  drawn  through  these  two  points  on  Figure  5-1 

It  will  be  noted  from  Figure  5-1  that  the  binomial 
distribution  points  lie  off  the  normal  distribution  line.   It 
can  therefore  be  seen  that  for  small  values  of  p  (in  the  range 
of  15%)  and  for  relatively  small  values  of  n,  the  use  of  the 
normal  distribution  would  not  be  suitable  to  simulate  the  bi- 
nomial distribution  for  the  entire  range  of  x^. 

With  this  in  mind  the  binomial  distribution  for  this 
study  can  be  simulated  indirectly  by  assigning,  for  each  value 
of  x,  a  series  of  consecutive  numbers  equal  to  the  values  of 
w(x)  listed  on  Table  5-1. 

TABLE  5-3 
SIMULATION  OF  BINOMIAL  DISTRIBUTION  -  RANDOM  NUMBERS 
Number  of  Successes  (x)    w(x)  x  1000   Random  Number  Ranges 

0  0388  0000  -  0387 

1  1368  0388  -  1755 

2  2293  1756  -  4048 

3  2428  4049  -  6476 

4  1821  6477  -  8297 

5  1028  8298  -  9325 

6  0454  9326  -  9779 

7  0160  9780  -  9939 

8  0046  9940  -  9985 

9  0011  9986  -  9996 
10  0003  9997  -  9999 
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If  we  now  generate  a  four  digit  random  number  from  some  random 
number  generator  this  number  can  be  equated  to  a  specific  value 
of  x.   See  Table  5-3.   As  an  example,  if  a  random  number  of 
9395  were  generated  this  would  be  equivalent  to  6  successes. 
Likewise,  if  a  number  0035  were  generated  this  would  equate  to 
0  successes,  or  all  dry  holes.   With  enough  trials  the  binomial 
distribution  can  be  simulated  in  this  manner. 

5 . 4  Coefficient  of  Variation 

From  Equations  6  and  7  the  coefficient  of  variation  can 
be  derived  as  follows: 

cv  =    Q     =  /  npq  (9) 

x0      np 
Since  q  =  (1-p) ,  Equation  9  can  be  written: 


fV  =  /  nP(x-p) 
uv       np 


From  which  we  obtain 


cv  =  JL  \/ii£  nrn 

/n~  V  P  U°) 

Using  the  data  for  number  of  wells  drilled  and  prob- 
ability of  success  previously  noted,  it  is  possible  to  calculate 
the  Coefficient  of  Variation  using  Equation  10  as  follows: 


CV  =  JL  l/^tP  =  "At?  (v/ST^D 
^20  V  .15     /2U 

CV  =  0.53  (11) 


Figure  5-1 
The  Normal  §  Binomial  Distributions 
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CHAPTER  VI 
OIL  RESERVOIR  SIZE 

It  goes  without  saying  that  the  economic  success  or 
failure  of  a  drilling  venture  may  well  be  determined  by  the 
size  of  the  field  discovered  as  a  result  of  a  successful  well. 
Of  course  the  size  of  the  field  is  unknown  prior  to  drilling. 
It  is  therefore  essential,  if  we  are  going  to  attempt  to  es- 
tablish the  relative  profitability  of  the  venture,  to  determine 
two  items.   First,  whether  there  is  a  general  mathematical  ex- 
pression which  can  be  used  to  characterize  the  probability 
distribution  of  field  sizes  in  the  geological  province.   In 
other  words,  do  the  field  sizes  vary  in  such  a  way  that  their 
probability  of  occurrence  can  be  predicted?   Secondly,  if  a 
general  field  size  distribution  function  can  be  established  as 
a  model,  what  are  the  boundary  conditions  we  should  establish 
for  our  particular  venture?   Each  of  these  questions  will  now 
be  explored. 

6 . 1  Field  Size  Distribution 

At  first  glance  it  would  seem  quite  unlikely  that  a 
general  mathematical  model  of  field  size  distributions  could  be 
established  -  considering  the  complex  process  of  oil  reservoir 
formation,  or,  if  by  chance  a  model  could  be  established  for 
one  geological  province  it  would  have  general  application  in 
another  province. 

Some  of  the  earliest  work  in  this  area  was  reported  by 
J.J.  Arps  and  T.G.  Roberts  [3]  in  1958  after  a  study  of  reserve 
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data  from  the  Denver- Julesburg  Basin.   In  their  study,  a  5.7 
million  acre  sample  was  chosen  on  the  east  flank  of  the  Basin 
(located  in  Nebraska,  Wyoming  and  Colorado)  and  the  reserve 
estimates  of  all  338  fields  found  to  that  date  in  the  Lower 
Cretaceous  sands  were  investigated.   Arps  and  Roberts  arranged 
the  fields  in  nine  groups  according  to  their  areal  extent  and 
estimated  the  ultimate  primary  recovery  from  each  group.   The 
limits  of  these  groups  were  chosen  in  such  a  way  that  each 
group  represented  fields  about  twice  as  large  as  the  preceeding 
one.   The  nine  groups  thus  arranged  are  represented  on  Table  6-1 

TABLE  6-1 
OIL  RECOVERY  DATA  ON  DENVER- JULESBURG  BASIN  FIELDS 


Range  of  Productive   Number  of   Estimated  Average 
Group         Acres  Fields    Ultimate  Recovery 

(M  bbls) 


1 

28 

2 

28  - 

57 

3 

57  - 

113 

4 

113  - 

226 

5 

226  - 

453 

6 

453  - 

905 

7 

905  - 

1810 

8 

1810  - 

3620 

9 

3620 

56 
77 
41 
70 
40 
33 
19 
1 
1 


10.1 

53.0 

165.0 

295.7 

957.6 

1,876.2 

4,481.6 

11,460.0 

50,750.0 


When  a  frequency  -  density  distribution  of  the  data 
presented  in  Table  6-1  was  plotted  versus  Average  Ultimate 
Recovery  on  logarithmic-probability  paper  a  straight  line  can 
be  best  fit  through  the  resulting  points.  (See  Figure  6-1). 
This  demonstrates  that  the  field  sizes  in  the  Denver- Julesburg 
Basin  are  Lognormally  distributed. 
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While  a  Lognormal  distribution  best  fits  the  data  in 
the  Denver-Julesburg  Basin  there  is  still  no  reason  to  general- 
ize and  say  that  all  reserves  follow  a  Lognormal  frequency  dis- 
tribution.  If,  however,  the  Lognormal  distribution  remains 
invariant  under  a  wide  variation  in  the  definition  of  reserves 
and  geological  regions  then  a  generalization  can  be  made. 

To  test  this  invariancy  Gordon  M.  Kaufman  [27]  ,  in  an 
extensive  study  reported  in  1962,  accumulated  data  from  several 
different  sources  including:   The  Oil  and  Gas  Conservation 
Board  of  Alberta  Report  to  the  Lieutenant-Governor- in-Council 
as  summarized  in  The  Financial  Post;  Oil  and  Gas  Journal  stat- 
istics on  two  states  for  1946  and  for  1960;  an  individual 
company's  estimate  of  reserves  in  four  basins;  and  the  Arps  - 
Roberts  article  on  the  Denver-Julesburg  Basin. 

The  data  accumulated  by  Kaufman,  covering  fields  in 

four  coastal  basins,  fields  in  North  and  South  Louisiana  (both 

in  1946  and  in  1959),  fields  in  Oklahoma,  etc.  were  plotted  on 

sixteen  (16)  separate  graphs  on  logarithmic-probability  paper. 

Kaufman,  in  analyzing  the  results  of  these  plots,  writes  [28] 

The  graphs  displaying  this  data  tentatively  suggest 
that: 

1.  Within  a  given  area  the  functional  form  of  empirical 
histograms  of  reported  field  sizes  does  not  change 

(a)  over  time, 

(b)  with  changes  in  the  definition  of  "reserves", 

(c)  with  changes  in  the  minimum  size  of  field 
reported,  or 

(d)  from  one  geographic  area  to  another. 

2.  The*  functional  form  is  the  same  for  a  well-defined 
geological  basin  such  as  the  Denver-Julesburg  as  for 
an  arbitary  geographic  area  such  as  a  state. 
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Although  each  of  /the  sixteen  (16)  graphs^/  is  a  plot 
of  data  from  fielcTs  of  varying  geological  ages  and  of 
varying  length  of  production  lives,  the  invariance  in 
the  functional  forms  is  striking  enough  to  warrant  the 
conjecture  that  added  data  blocked  on  these  two  factors 
would  not  alter  the  conclusions  stated  above. 

A  visual  observation  of  these  graphs  reveals  that  most  points 

generally  fall  along  a  straight  line  indicating  a  Lognormal 

frequency  distribution  of  field  size. 

In  another  analysis  of  oil  field  size  distributions, 
Folkert  Brons  reports  the  result  of  a  study  of  the  predicted 
ultimate  recovery  of  75  Southern  Louisiana  Miocene  Fields  [10]. 
As  shown  on  Figure  6-2,  when  the  data  is  plotted  on  logarithmic 
probability  paper  the  points  fall  along  a  straight  line  indi- 
cating a  Lognormal  distribution  function  for  these  fields. 

On  the  basis  of  the  evidence  presented  in  the  studies 
cited  above,  of  reservoirs  with  varying  reserve  sizes,  geologic 
origin  and  physical  characteristics,  it  will  be  assumed  for 
further  analysis  that  a  Lognormal  distribution  function  can 
serve  as  a  general  mathematical  model  for  reservoir  size 
distributions . 

6. 2  The  Lognormal  Distribution 

The  lognormal  distribution  can  be  defined  as  the  dis- 
tribution of  a  random  veriable  whose  distribution  obeys  a 
normal  distribution  function.   Thus,  if  x  is  lognormally  dis- 
tributed then  y  =  In  x  is  normally  distributed  and, 

f-(y-yo)2" 


w(y)  =  — -— r  exp 

0"y  /2tt 


V 


38 


39 


If  the  distribution  is  desired  in  terms  of  x  rather 

than  y  a  change  in  variables  is  performed  which  results  in  the 

probability  density  function  for  y  expressed  in  terms  of  x  as 

follows : 

2^ 


w(x)  = 


ay   /2? 


exp 


-(In  x-yQ) 


2a 


y 


(2) 


where:  x  =  independent  variable 
ay  =  standard  deviation 
y0  =  mean  value  of  y 
The  equation  gives  a  lognormal  distribution  curve, 
with  its  characteristic  positive  skew,  as  shown  below. 


Lognormal  Distribution  -  Positive  Skew 
Figure  6-3 
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Since  y  is  normally  distributed: 

y  (median  =  y  (mean)  =  y(mode),  (3) 

but  since  x  =  e^,  x  is  not  normally  distributed  and  the  values 
of  x  (median) ,  x  (mean)  and  x(mode)  must  be  calculated  separ- 
ately [9]  : 

The  Median,  xm: 

y  =  In  x 

xm  =  Sm  (4) 

m 

since  ym  =  yQ 

xm  =  e  (4A) 


and,  In  xm  =  yQ  (4B) 


The  Mean,  xQ : 

2 
y   +  Oy'' 

xn  =  e  °    ~T  (5) 


■o 


or,  In  xn  =  y, 


+ 


a 


X 


2 


ko    7o     2  (5A) 

The  Mode,  xmode: 


^o  -   2 


mode  =  e    °y  (6) 


or>  ln  xmode  =  y0    ~   °y2  (6A) 
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2 
Further,  the  Variance ,  a        is  given  by  the  following  expression 


=  exp  {  2yQ  +  a   2    }  (e^  -1) 


x 


(7) 


Inserting  Equation  5  in  Equation  7  we  obtain: 

a   2  =  x„2  (e°y2-l)  (8) 


From  which  the  Coefficient  of  Variation,  CV,  may  be  expressed 

as : 

(9) 

CV  =  —  =   eay  -1 
xo 

For  small  values  of  ay  we  can  approximate  Equation  9  as: 

CV  =  ay  (10) 


6 . 3  Boundry  Conditions 

If  we  assume  that  the  field  (or  reservoir)  sizes  in 
any  area  are  lognormally  distributed  and  will  fall  along  a 
straight  line  when  plotted  on  logarithmic-probability  paper 
then  the  next  problem  is  one  of  specifing  the  position  of  the 
line  which  would  represent  the  reservoir  size  distribution  in 
the  specific  location  where  the  drilling  venture  is  to  be 
conducted. 

This  line  could  be  drawn  if  someone  could  predict  the 
median  value  of  the  reserve  size  to  be  discovered  by  a  (yet 
undrilled)  successful  well  (i.e.  the  reserve  size  corresponding 
to  a  50%  probability)  and  the  standard  deviation  of  the  dis- 
tribution representing  the  logarithm  of  the  field  sizes. 
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Then,  by  inserting  these  values  in  Equations  (4B) ,  (5A)  and 
(6A)  previously  noted  for  median,  mean  and  mode  a  straight  line 
could  be  drawn  on  log-probability  paper.   The  problem  is,  how- 
ever, that  the  operator  (or  his  agent)  would  probably  not  be 
able  to  specify  either  the  median  or  the  standard  deviation. 
What  he  most  likely  would  be  able  to  specify  is  his  estimate  of 
the  highest  and  lowest  values  of  reserves  he  would  expect  with, 
say,  90%  certainty.   He  would  probably  base  these  estimates  on 
either  his  own  experience  or  the  experience  of  others  who  have 
drilled  in  the  area.   As  an  example,  the  operator  may  not  be 
able  to  specify  an  exact  reserve  size,  although  he  might  say 
that  he  would  expect  the  reserves  to  run  from  5  million  to  150 
million  barrels  with  90%  confidence.   If  these  two  points  were 
plotted  on  log-probability  paper  -  with  the  5  million  barrel 
figure  plotted  at  the  5%  probability  point  and  the  150  million 
barrel  figure  plotted  at  the  951  probability  point  a  straight 
line  could  be  drawn  between  them  which  would  represent  the 
Lognormal  distribution  of  all  reserves  in  the  area. 

With  the  line  plotted,  the  median  value  of  Q  (or  xm  in 
Equation  4B) ,  corresponding  to  a  501  probability,  and  the  value 
of  Q  median  plus  one  standard  deviation,  corresponding  to  an 
84%  probability,  can  be  read  directly.   In  addition,  to  use  this 
model  in  our  determination  of  the  relative  profitability  of  the 
venture  using  the  Monte  Carlo  technique,  it  is  essential  that 
we  be  able  to  determine  a  reservoir  size  and  corresponding  value 
of  probability  (that  the  reservoir  size  will  not  be  exceeded) 
for  each  point  along  the  line. 


43 


To  illustrate  this  procedure,  as  well  as  to  generate 
data  for  use  in  the  Monte  Carlo  studies  conducted  as  part  of 
this  thesis,  let  us  assume  we  are  going  to  drill  for  oil  on- 
shore in  the  Southern  Louisiana  area.   By  going  to  the  approp- 
riate issue  of  the  Oil  and  Gas  Journal  [48]  we  would  find  a 
listing  of  the  sizes  of  all  fields  in  the   Southern  Louisiana 
on-shore  area.   These  fields  could  then  be  arranged  in  ascend- 
ing order  of  size  (Q)  and  a  relative  frequency  density,  w,  and 
cumulative  frequency,  W(Q  £  Q.)  calculated.   This  tabular 
listing,  along  with  the  calculation  of  relative  and  cumulative 
frequency  is  shown  as  ACTUAL  Values  in  Table  6-3. 

Having  the  actual  values  of  cumulative  frequency  cal- 
culated, these  values  can  be  plotted  on  log-probability  paper 
with  the  ordinate  being  the  cumulative  frequency  and  the  ab- 
scissa being  the  field  size.   The  triangular  points  on  Figure 
6-4  represent  a  plot  of  the  actual  data  for  the  Southern  Louis- 
iana on-shore  area. 

The  operator  can  now  use  these  points  (or  any  other 
data  he  has  available)  as  a  guide  in  drawing  a  straight  line 
on  the  log-probability  paper  beginning  at  some  upper  value  of 
Q,  corresponding  to  a  95%  probability,  and  at  some  lower  value 
of  Q,  corresponding  to  a  5%   probability.   In  this  example  the 
points  chosen  were  79  million  barrels  as  the  upper  value  and 
5  million  barrels  as  the  lower  value.   With  this  line  drawn 
(Curve  A,  Figure  6-4)  it  can  now  be  said  that,  for  the  purposes 
of  this  profitability  analysis,  90%  of  all  reserves  discovered 
by  a  successful  well  will  fall  along  the  straight  line,  5%  of 
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them  will  exceed  the  upper  value  and  5%   will  be  less  than  the 
lower  value.   Theoretically  then,  there  is  a  5%   possibility 
that  a  reservoir  could  be  discovered  with  infinite  reserves  on 
one  extreme  and  a  5%  possibility  of  a  reservoir  with  a  0  res- 
erve on  the  other.   Practically  of  course,  this  is  not  reason- 
able since  0  reserve  means  no  oil  and  an  unsuccessful  well. 
If  this  straight  line  represents  the  reserve  size  distribution 
for  a  successful  well  then  some  lower  limit  must  be  established 
which  would  represent  the  reserve  size  below  which  the  well  in 
question  would  not  be  considered  successful,  even  though  some 
oil  was  discovered.   In  this  example  a  lower  limit  of  1  million 
barrels  was  chosen.   Therefore,  we  say  to  be  considered  succ- 
essful a  well  must  discover  at  least  a  1  million  barrel  reserve 

At  the  upper  extremes  it  would  be  foolish  to  say  that 
an  infinite  reservoir  could  be  discovered  even  though  this  is 
possible  by  an  extrapolation  of  Curve  A,  Figure  6-4,  to  a  1001 
probability  point.   Realizing  that  this  is  not  physically 
possible  an  arbitrary  upper  limit  should  be  established.   In 
this  example  it  was  observed  that  in  the  Southern  Louisiana 
on-shore  area  the  maximum  field  size  is  137  million  barrels 
(see  Table  6-3,  ACTUAL).   On  this  basis  an  upper  limit  of  140 
million  barrels  was  chosen. 

6 . 4  The  Lognormal  Distribution  of  Field  Sizes 

With  the  boundary  conditions  now  established  and  Curve 
A, Figure  6-4  drawn,  it  is  now  possible  to  proceed  with  the 
calculation  of  the  mean,  median  and  mode  as  follows: 


The  Median,  x  : 


Q  at  50%  =  20.0  million  barrels 


let  Median  =  xi 
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From  Equation  4B 


=  y„  =  In  x 


m 

y0 


m 
In  20.00 

2.996 


y0  +  a   =  In  x  (at  W  <  x  = 

y0  +  oy  =  In  46.0 

2. 996  +  o  =  3.829 

a  =  0.833 


84%) 


oyc    =  0.694 


a  2  =  0.347 


The  Mean,  x, 


From  Equation  5A 

In  x( 
In  x 


+  °-X 


2 


^0  '  T 

2.996  +  0.347 


In  x„  = 


3.343 
28.3 


million  barrels 
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The  Mode,  xmode: 


From  Equation  6A 


In  x„,rt-,   =  y_  -  a,,2 
mode   7o    y 


ln  xmode  =  2-996  "  °'694 


lnxm  ,   =  2.302 
mode 


In  xmQ(je  =  10.0     million  barrels 


From  Equation  9 


.a...  ^7T 


cv 

xo 


CV   =     \/e°-694-l 


CV  =  /2. 0017-1  =  /1.0017 

CV  =  1.00 

If  we  assign  random  numbers  to  represent  W(Q  <,  Q^)  we 
can,  from  Curve  A,  Figure  6-4,  manually  determine  the  corres 
ponding  value  of  Q  as  follows: 
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TABLE  6-2 
FIELD  SIZES  -  RANDOM  NUMBERS 

W(Q  <  Qi)  Q 

Random  No.  (in  %)  (x  106bbl) 


05  05  5.0 

06  06  5.5 

07  07  5.8 


93  93  69.0 

94  94  73.0 

95  95  80.0 


While  the  process  of  manually  reading  the  graph,  tab- 
ulating the  values  of  Q  corresponding  to  Random  Variables 
which  represent  W(Q  £  Q-)  and  inputing  these  numbers  as  data 
in  the  Monte  Carlo  computer  program  is  possible,  the  process 
is  quite  laborious.   To  equip  ourselves  for  an  efficient  sim- 
ulation system  a  Lognormally  Distributed  Variable  Generator 
can  be  used  to  generate  the  values  of  Q  on  a  high  speed 
digital  computer. 

Such  a  Generator  was  initially  developed  by  McMillan 
and  Gonzalez  [34].   A  version  of  the  Generator,  modified  for 
use  in  this  problem,  incorporating  the  concepts  of  upper  and 
lower  limits,  is  included  in  Appendix  A.   Using  this  Generator 
it  is  necessary  only  to  input  the  median  Value  of  Q  and  the 
value  of  Q  at  W(Q  £  Qi)  =  841,  both  easily  obtainable  from 
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the  plot  of  a  straight  line  connecting  the  5%  and  95%  points, 
and  the  upper  and  lower  limits,  as  previously  discussed. 
With  this  data  inputed,  the  Generator  will  provide  field  sizes 
(Q)  which  are  lognormally  distributed  within  the  boundary 
conditions  established  for  the  program. 

In  order  to  test  the  reliability  of  the  deviates 
provided  by  the  Lognormal  Generator,  1000  values  of  Q  were 
generated  using  the  program  included  in  Appendix  A  and  the 
boundary  conditions  previously  discussed.   The  results  of  this 
test,  along  with  a  calculation  of  relative  and  cumulative 
frequency,  is  included  on  Table  6-3  under  the  heading  of 
SIMULATION.   For  comparison  purposes  the  values  provided  by 
the  Generator  were  placed  beside  the  corresponding  actual 
values.   Additionally,  Curve  B,  Figure  6-4  was  plotted  with 
points  generated  by  the  Lognormal  Generator. 

By  comparing  Curves  A  and  B  and  the  actual  data  points 
on  Figure  6-4  it  will  be  seen  that  excellent  correlation 
exists  except  in  the  lower  values  of  Q  -  but  even  there  the 
difference  is  not  considered  significant.   On  this  basis, the 
Lognormally  Distributed  Variable  Generator  will  be  used  to 
simulate  lognormally  distributed  values  of  field  size  in  this 
thesis . 


Figure  6-4 
Cumulative  Distribution  of  Field  Sizes 
Southern  Louisiana  On-Shore 


Field  Size  -  Q  (xl06bbl) 
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TABLE  6-3 


SIZE  OF  FIELD  -  SOUTHERN  LOUISIANA  ON-SHORE 
(1970  Estimates) 


ACTUAL 

SIMULATION 

Field 

Freq. 

Rel . Freq . 

Cum. Freq. 

Freq. 

Rel . Freq . 

Cum. Freq . 

Size-Q 

£ 
(w=n-xl00) 

f 
(w=ttx100) 

(xl06bbl) 

(f) 

(W  <  Q) 

(f) 

(W  <  Q) 

m 

m 

w 

m 

1 

3 

0.3 

0.3 

2 

11 

1.1 

1.4 

3 

18 

1.8 

3.2 

4 

25 

2.5 

5.7 

5 

24 

2.4 

8.1 

6 

2 

2.4 

2.4 

33 

3.3 

11.4 

7 

3 

3.5 

5.9 

27 

2.7 

14.1 

8 

2 

2.4 

8.3 

34 

3.4 

17.5 

9 

0 

8.3 

19 

1.9 

19.4 

10 

5 

5.9 

14.2 

45 

4.5 

23.9 

11 

4 

4.7 

18.9 

34 

3.4 

27.3 

12 

5 

5.9 

24.8 

25 

2.5 

29.8 

13 

0 

24.8 

23 

2.3 

32.1 

14 

5 

5.9 

30.7 

37 

3.7 

35.8 

15 

7 

8.2 

38.9 

28 

2.8 

38.6 

16 

3 

3.5 

42.4 

37 

3.7 

42.3 

17 

4 

4.7 

47.1 

32 

3.2 

45.5 

18 

5 

5.9 

53.0 

25 

2.5 

48.0 

19 

2 

2.4 

55.4 

25 

2.5 

50.5 

20 

1 

1.2 

56.6 

27 

2.7 

53.2 

21 

1 

1.2 

57.8 

21 

2.1 

55.3 

22 

1 

1.2 

59.0 

18 

1.8 

57.1 

23 

4 

4.7 

63.7 

22 

2.2 

59.3 

24 

1 

1.2 

64.9 

20 

2.0 

61.3 

25 

0 

64.9 

11 

1.1 

62.4 

26 

1 

1.2 

66.1 

19 

1.9 

64.3 

27 

1 

1.2 

67.3 

15 

1.5 

65.8 

28 

3 

3.5 

70.8 

12 

1.2 

67.0 

29 

1 

1.2 

72.0 

13 

1.3 

68.3 

30 

1 

1.2 

73.2 

20 

2.0 

70.3 

31 

1 

1.2 

74.4 

15 

1.5 

71.8 

32 

74.4 

8 

0.8 

72.6 

33 

74.4 

13 

1.3 

73.9 

34 

74.4 

18 

1.8 

75.7 

35 

2 

2.4 

76.8 

6 

0.6 

76.3 

36 

1 

1.2 

78.0 

13 

1.3 

77.6 

37 

78.0 

4 

0.4 

78.0 

38 

78.0 

8 

0.8 

78.8 

39 

78.0 

7 

0.7 

79.5 

40 

2 

2.4 

80.4 

6 

0.6 

80.1 

TABLE    6-3      Continued 


51 


ACTUAL 

SIMULATION 

Field 

Freq. 

RelFreq . 

C  um . F  r  e  q . 

Freq . 

Rel . Freq . 

Cum. Freq . 

Size-Q 

£ 
(w=irxl00) 

f 

(xl06bbl) 

(f) 

(W  <  Q) 

(f) 

(w=-xl00) 

(W  <  Q) 

(%) 

(T) 

m 

m 

41 

1 

1.2 

81.6 

5 

0.5 

80.6 

42 

2 

2.4 

84.0 

8 

0.8 

81.4 

43 

84.0 

9 

0.9 

82.3 

44 

84.0 

3 

0.3 

82.6 

45 

1 

1.2 

85.2 

6 

0.6 

83.2 

46 

85.2 

8 

0.8 

84.0 

47 

1 

1.2 

86.4 

4 

0.4 

84.4 

48 

86.4 

9 

0.9 

85.3 

49 

86.4 

4 

0.4 

85.7 

50 

86.4 

2 

0.2 

85.9 

51 

1 

1.2 

87.6 

9 

0.9 

86.8 

52 

87.6 

5 

0.5 

87.3 

53 

1 

1.2 

88.8 

5 

0.5 

87.8 

54 

88.8 

2 

0.2 

88.0 

55 

88.8 

6 

0.6 

88.6 

56 

88.8 

5 

0.5 

89.1 

57 

88.8 

2 

0.2 

89.3 

58 

88.8 

7 

0.7 

90.0 

59 

1 

1.2 

90.0 

3 

0.3 

90.3 

60 

2 

2.4 

92.4 

6 

0.6 

90.9 

61 

92.4 

6 

0.6 

91.5 

62 

92.4 

2 

0.2 

91.7 

63 

92.4 

7 

0.7 

92.4 

64 

1 

1.2 

93.6 

2 

0.2 

92.6 

65 

93.6 

1 

0.1 

92.7 

66 

93.6 

5 

0.5 

93.2 

67 

93.6 

1 

0.1 

93.3 

68 

93.6 

2 

0.2 

93.5 

69 

93.6 

4 

0.4 

93.9 

70 

93.6 

3 

0.3 

94.2 

71 

93.6 

3 

0.3 

94.5 

72 

93.6 

3 

0.3 

94.8 

73 

93.6 

94.8 

74 

93.6 

3 

0.3 

95.1 

75 

1 

1.2 

94.8 

1 

0.1 

95.2 

76 

94.8 

2 

0.2 

95.4 

77 

94.8 

1 

0.1 

95.5 

78 

94.8 

1 

0.1 

95.6 

79 

94.8 

1 

0.1 

95.7 

80 

94.8 

95.7 

81 

94.8 

1 

0.1 

95.8 

82 

94.8 

3 

0.3 

96.1 

TABLE    6-3      Continued 
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ACTUAL 

SIMULATION 

Field 

Freq . 

Rel . Freq . 

Cum. Freq . 

Freq. 

Rel . Freq . 

Cum. Freq. 

Size-Q 

(w=^xl00) 

(w=£xl00) 

(xl06bbl) 

(f) 

(W  <  Q) 

(f) 

(W  <  Q) 

m 

m 

w 

(%) 

83 

94.8 

96.1 

84 

94.8 

1 

0.1 

96.2 

85 

94.8 

4 

0.4 

96.6 

86 

94.8 

96.6 

87 

94.8 

96.6 

88 

94.8 

2 

0.2 

96.8 

89 

94.8 

96.8 

90 

94.8 

1 

0.1 

96.9 

91 

94.8 

96.9 

92 

94.8 

96.9 

93 

1 

1.2 

96.0 

1 

0.1 

97.0 

94 

96.0 

3 

0.3 

97.3 

95 

1 

1.2 

97.2 

4 

0.4 

97.7 

96 

97.2 

2 

0.2 

97.9 

97 

97.2 

97.9 

98 

97.2 

97.9 

99 

97.2 

97.9 

100 

97.2 

97.9 

101 

97.2 

2 

0.2 

98.1 

102 

97.2 

1 

0.1 

98.2 

103 

97.2 

1 

0.1 

98.3 

104 

97.2 

98.3 

105 

97.2 

1 

0.1 

98.4 

106 

97.2 

1 

0.1 

98.5 

107 

97.2 

98.5 

108 

1 

1.2 

98.4 

98.5 

109 

98.4 

98.5 

110 

1 

1.2 

99.6 

98.5 

111 

99.6 

98.5 

112 

99.6 

98.5 

113 

99.6 

1 

0.1 

98.6 

114 

99.6 

1 

0.1 

98.7 

115 

99.6 

98.7 

116 

99.6 

98.7 

117 

99.6 

1 

0.1 

98.8 

118 

99.6 

98.8 

119 

99.6 

1 

0.1 

98.9 

120 

99.6 

98.9 

121 

99.6 

98.9 

122 

99.6 

98.9 

123 

99.6 

3 

0.3 

99.2 

124 

99.6 

99.2 

125 

99.6 

1 

0.1 

99.3 
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ACTUAL 

SIMULATION 

Field 
Size-Q 

(xl06bbl) 

Freq.   Rel.Freq.   Cum.Freq. 
(£)    (w=|ocl00)    (W  <  Q) 

Freq.   Rel.Freq.   Cum.Freq. 

(£)    (w=|xl00)    (W  <  Q) 
(%)        (%) 

126 

127 
128 
129 
130 
131 
132 
133 
134 
135 
136 
137 
138 
139 
140 

99.6 
99.6 
99.6 
99.6 
99.6 
99.6 
99.6 
99.6 
99.6 
99.6 
99.6 
1       1.2      100.0 

100.0 
100.0 
100.0 

99.3 
99.3 
99.3 
99.3 
99.3 
99.3 
99.3 
99.3 

1        0.1        99.4 

99.4 
99.4 
99.4 

1       0.1       99.5 

99.5 

5       0.5      100.0 
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CHAPTER  VII 
DISCOUNTED  PROFIT  PER  BARREL 

7 . 1  General 

In  a  very  real  sense,  the  ultimate  profitability  of  an 
exploratory  well  drilling  venture,  where  oil  is  discovered,  is 
dependent  not  only  on  the  size  of  the  new  reserves  but  also 
(and  perhaps  more  importantly)  on  how  much  profit  the  operator 
will  receive  for  each  barrel  of  oil  produced.   The  case  can  be 
visualized  where  the  operator  in  a  successful  exploratory  drill 
ing  venture  discovers  vast  potential  oil  reserves  but  the 
expence  of  production  far  exceeds  the  posted  price  he  can 
receive  for  the  oil  produced.   In  this  case  no  matter  how  large 
the  potential  reserve,  the  profitability  of  the  venture  is  less 
than  0,  or  if  produced,  the  operator  would  suffer  a  loss.   The 
cases  of  off-shore  drilling  in  deep  water  and  production  on  the 
North  Slope  of  Alaska  (the  latter  receiving  wide  publicity) 
illustrate  this  notion.   It  is  therefore  essential  in  a  profit- 
ability analysis  to  estimate  the  expected  profit  per  barrel  on 
reserves  discovered  by  a  successful  well.   In  practice,  though, 
it  is  rare  that  the  operator  would  commit  himself  to  "the" 
expected  profit  per  barrel  for  the  venture.   He  could  perhaps 
provide  his  estimates  of  the  lowest,  highest  and  most  likely 
values  of  profit  per  barrel  basing  these  estimates  on  his 
experience  in  the  area  or  on  factors  similar  to  those  which 
permitted  him  to  estimate  the  probability  of  success  (see 

Chapter  V) . 
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7 . 2  Discounted  Profit  Per  Barrel  Estimate 

In  estimating  the  highest  lowest  and  most  likely  values 
of  profit  per  barrel  the  operator  would  first  have  to  provide 
his  estimate  of  the  highest,  lowest  and  most  likely  posted 
price  for  the  oil.   From  this  he  would  have  to  deduct  his  ex- 
penses including,  royalties,  severance  taxes,  operating  costs, 
drilling  and  completion  costs,  etc.   This  would  then  provide 
his  net  profit  per  barrel  for  each  or  the  three  estimates. 
But,  since  reserves  cannot  be  produced  all  at  once,  the  profit 
per  barrel  can  not  be  figured  on  the  basis  of  current  sales 
prices.   Production  will  be  spread  into  the  future  (say  10  years) 
and  future  money  is  not  as  valuable  as  present  day  money.   To 
account  for  this  financial  fact-of-life  the  operator  will  have 
to  convert  his  net  profit  per  barrel  calculated  at  today's 
prices  into  a  discounted  cash  flow.   In  this  discounting  (or 
deferment)  process  a  discounting  factor  is  developed  by  calcu- 
lating the  Present  Worth  of  a  $1  annuity  (i.e.  constant  income 
or  expenditure)  received  each  year  for  a  specified  length  of 
time,  n,  at  a  specified  discounted  rate,  i,  [16].  This  is  given 
by  the  expression: 


P  =  R 


["(1  +  i)n  -1 


(1) 


where:   P  is  the  Present  Worth 

R  is  the  amount  of  the  annuity  (taken  as 
$1  in  this  development) 
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n 


(1  +  i)"  -1 
1(1  +  i) 


n 


is  the  uniform  series  Present  Worth 
or  Present  Value  Factor. 


1  - 


The  Present  Worth  Factor  can  also  be  written  as 


or  Values  of  the  Present  Worth  can  be 


(1  +  i)n 


obtained  from  Tables.  If  these  values  are  then  divided  by  the 
number  of  discounting  periods,  n,  the  deferment  or  discounting 
factor  is  obtained. 

To  illustrate,  the  discount  factor  based  on  a  discount 
rate  of  101  for  a  period  of  10  years  is  calculated  as  follows: 

From  Equation  1, 
P  =  R 


(1  +  i) 


:  ^n 


i(l  +  i)n 

Let  R  =  $1 
i  =  10% 
n  =  10  years 
From  Tables  [17]  the  Present  Worth  Factor  is  6.1444 

P  =  1(6.1444) 
P  =  6.1444 

(2) 


p 
Discount  factor  = 


n 


6  144 
Discount  factor  =  z~= — 

Discount  factor  =  .6144 


Concerning  this  technique,  Grayson  [19]  writes: 

Many  operators,  ...  object  to  the  present  value  concept 
First,  there  are  those  who  scoff  at  the  idea  of  making 
any  present  value  computations  on  exploratory  ventures. 
There  are  already  so  many  uncertainties,  is  their 
belief,  that  to  refine  the  calculation  any  further  is 
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like  tilting  at  windmills.   The  answer  to  this  has 
been  stated  several  times:   a  present  value  calculation 
is  an  estimate,  but  an  educated  one  that  explicitly 
recognizes  time  as  a  realistic  decision  consideration. 
To  ignore  the  effect  of  time  or  to  try  to  cope  with  it 
in  an  implict  way  is  to  incur  the  danger  that  it  may 
be  overlooked  entirely  or  its  effect  miscalculated. 

On  the  basis  of  the  concepts  discussed  in  this  section, 
the  highest,  lowest  and  most  likely  values  of  profit  per  barrel 
can  be  computed.   Table  7-1  illustrates  this  process  and  is 
based  on  the  hypothetical  expenses  and  prices  for  a  drilling 
venture  in  the  Southern  Louisiana  on-shore  area. 

The  values  developed  in  Table  7-1  will  be  used  as  the 
base  figures  for  the  sensitivity  studies  contained  in  subsequent 
Chapters . 

7 . 3  The  Triangular  and  Beta  Distributions 

It  was  seen  in  Section  7.2  that  while  the  operator 
would  probably  be  unable  to  provide  the  exact  profit  per  barrel 
of  oil  produced  he  could  provide  his  estimates  of  the  highest 
profit  per  barrel  he  would  expect,  the  lowest  profit  per  barrel 
and  the  profit  per  barrel  he  feels  most  likely  to  occur.   It 
is  assumed  of  course  that  these  figures  will  be  discounted  to 
reflect  the  time  value  of  money.   Given  this  data,  the  question 
then  becomes,  what  distribution  function  best  suits  this  set  of 
parameters?   If  all  that  is  known  is  the  upper  limit,  the  lower 
limit  and  most  likely  value  (mode)  it  is  possible  to  describe 
this  knowledge  by  means  of  either  the  triangular  distribution 
function  or  the  beta  distribution  function.   Of  the  two,  the 
triangular  function  is  the  simplest  and  will  be  discussed  first. 


MOST 

LIKELY 

HIGH 

LOW 

$3.25 

3.50 

3.00 
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TABLE  7-1 
DISCOUNTED  PROFIT  PER  BARREL  -  SOUTHERN  LOUISIANA  ON-SHORE 

ITEM 

POSTED  PRICE 

DEDUCTIONS: 

ROYALTY: 

@l/6  .54  .50 

@l/8  .44 

SEVERANCE  TAXES  (see  Note  1)       .24      .26      .22 

OPERATING  COSTS  (including  lift 
cost,  labor,  lease  maint., 
etc.)  .10      .11     .09 

HOME  OFFICE  COSTS  OR  OVERHEAD 
(§  5%) 

GEOPHYSICAL  COSTS 

DRILLING  $  COMPLETION 
(see  Note  2) 

LEASE  ACQUISITION  (see  Note  3) 

TOTAL  DEDUCTIONS 

NET  PROFIT 

DISCOUNTED  PROFIT  (based  on 
discounting  at  10%  for  10 
years:   discount  factor  = 
.6144) 

NOTE  1:  [49]  a.   State  and  local  severance  taxes  paid  in  1968 

=  $233,017,000. 

b.  Total  Petroleum  liquid  production  in  1968  = 
2,526,008  bbl/day  (average). 

c.  Average  field  price  of  oil  (1969)  =  $3.33 


.16 

.18 

.15 

.01 

.01 

.01 

.40 

.40 

.40 

.02 

.02 

.02 

-$1.82 

-1 

.62 

-2.19 

1.43 

1 

.88 

.81 

$.88 

$1 

.15 

$.50 

$233,071,000 *  9C7/KK1 

2,526,008  bbl/day  x  365   "  *-^/&DJ- 
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TABLE  7-1  Continued 


however,  at  $3.30/bbl 

Severance  taxes  =  : ^52  =  7  351 

3.30 

say  7.4%  of  posted  price 


NOTE  2: [50]   a 
b 


Assume  production  at  150  bbl/day. 

Estimated  Cost  of  drilling  and  equipping  wells 
in  1970  =  $196,500.   Allowing  for  inflation 
§   deeper  wells  assume  cost  of  drilling  = 
$217,500/well. 

Assume  recovery  of  217,500  in  10  years  then 
recovery/year  =  $21,750. 

$  2 1  7^0 
150  bbl/day  x  365  =  $-40/bbl 


NOTE  3:   a. 


b. 


c. 


d. 


Assume  Lease  Costs  at  $60/acre  where  a  completion 
is  made. 

Assume  well  spacing  at  40  acres; 

then  Lease  Costs  $2400/well  where  well  is  completed 

Assume  lease  on  unproductive  land  to  be  equivalent 

to  $2400/well; 

then  total  lease  costs  =  $4,800/well. 

Assume  recovery  in  5  years  with  production  rates 
=  150/bbl/well. 

$4800   *Q,.n/   ,, 
— g —  =  $960/well 


$960 
150  bbl/day  x  365  "  'U1/ 


say  $.02/bbl 


60 


7.4  The  Triangular  Distribution  Function 

The  triangular  distribution  is  defined  by  three  values 
of  the  variable;  i.e.  the  lowest  valur  of  x  =  a,  the  highest 
value  of  x  ■  c  and  the  most  likely  value  of  x  =  b.   The  prob- 
ability distribution  function  w(x)  can  take  three  forms  as 
shown  below: 


w(x) 


w(b) 


■  x 


Figure  7-1A  Triangular  Distribution  -  Positive  Skew 
w(x) 


*-    X 


Figure    7-1B  Triangular  Distribution-    Symmetrical 
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w(b) 


Figure  7-1C  Triangular  Distribution  -  Negative  Skew 


Since   /  w(x)  dx  =  1  (=  the  area  of  the  triangle)  for  a 
a 

continuous  probability  distribution: 


/  w(x)  dx  =  l/2w(b) (c-a)  =  1 


from  which 


w(b)  = 


c-a 


(3) 


(4) 


In  unpublished  lecture  notes  Folkert  Brons  [11] 
proves  the  following  properties  for  the  triangular 
distribution: 


mean,  x 


a  +  b  +  c 


mode,  w(b) 


median,  x  =  c- 

'   m 


c-a 


(c-a)(c-b) 


(5) 
(6) 
(7) 


for  positive  skew  (see  Figure  7-1A) 
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xm  =  a  +  ^(c-a)(b-a)  (g) 


for  negative  skew  (see  Figure  7-lc) 
Brons  has  also  derived  the  cumulative  distribution  function 
[11]  and  has  demonstrated  that: 
For  x  ^  b 

(x^a)2 
W<*  *   *0    =  (b-aHc-a)  O) 

or,   xi  =  a  +  V  W(x  <  x^  (b-a)  (c-a)       (10) 

For  x  £  b 

wr   >   .      Cc-Xj)2  (11) 

W(X  >  X-)  =  -f r-TT1 r 

i     (c-b) (c-a) 
or,   xi  =  c  -  ^W(x  >  xi)  (c-b)  (c-a)        (12) 

Further, 

w(*  ib)  -  |^|  (13) 

The  variance  for  the  triangular  distribution  is  given  by: 

7     2     2 
2   a   +  b  +  c   -  ab  -   ac  -  be 
a  -  fg  (14) 

From  which  the  Coefficient  of  Variation  may  be  derived  as 
follows : 


2,22, 
a   +  b   +  c   -  ab  -ac  -  be 


rv    ax  .  J 18 

cv  _  — _n5__ 
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2        ,2  2  , 

wa      +   b      +   c      -    ab    -    ac    -    be 

CV   =  a   +    b   +    c 


CV   = 


_   {JTV+ 


c      -    ab    -    ac    -    be 


1.4l4(a   +   b   +   c)  (15) 


7 . 5  Simulation  of  the  Triangular  Distribution 

Using  the  data  for  the  highest,  lowest  and  most  likely 
values  of  discounted  profit  per  barrel  previously  developed  it 
is  possible  to  simulate  the  cumulative  triangular  distribution 
of  profit  per  barrel  W(x)  as  follows: 

Let  a  =  $.50 
b  =  $.88 
c  =  $1.15 

x  =  specific  values  of  expected  discounted 
profit  per  barrel 
then  for,  $ . 50  <  xi  <  $ . 88 


x-  =  a  + 


\/w(x  <  xi)(b-a)(c-a) 
:.  =  .50  +  \/w(x  <  x^  (.88-.50)(1.15  -.50) 
:i  =  .50  +  Vw(x  <  x^  (.2470)  (16) 


for,  $1.15  >  x.  >  $.88 
'       —   l  — 


x.  =  c  -  \/w(x  >  xi)  (c-b)  (c-a) 


:±   =  1.15  -  \/w(x  >  xi)(1.15-.88)  (1.15-.50) 


:.  =  1.15  -  Vw(x  >  xi)(.1755)  (17) 
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and, 


W(x  <  b)  = 


c-a 


.88-. 50 
w(x  4  $.88)  =  i.is-.so 

W(x  <  $.88)  =  .585  (18) 

Further,  from  Equation  15  the  Coefficient  of  Variation  may  be 
calculated  as  follows: 


ya  +  b   +  c   -  ab  -  ac  -  be 
LV  "     1.414  (a  +  b  +  c)  (15) 

rv    V(.50)2+(.88)2+(1.15)2-(.50)(1.15)-(.88)(1.15) 
cv  ~  1.414  (.50  +  .88  +  1.15) 


CV  =  -^51 
LV    3.580 

CV  =  .157 

A  series  of  random  numbers  can  now  be  assigned  to  represent 

values  of  W(x  <  xi)  from  .0000  to  .999999 From 

Equation  16,  if  this  random  number  was  .585  or  less  the  random 
number  would  be  substituted  into  Equation  10,  to  obtain  a 
value  of  x..   If  the  random  number  was  .585  or  greater 

W(x  <  xi)  =  1.0  -  W(x  >  xi)  (19) 

would  be  used  to  obtain  W(x  >  x^)    and  then  that  value  would 
be  substituted  into  Equation  12,  to  obtain  the  corresponding 
value  of  x^.   This  process  is  illustrated  in  Table  7-2. 
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TABLE  7-2 
ILLUSTRATION:   PROFIT  PER  BARREL  ASSUMING  TRIANGULAR 
DISTRIBUTION 


Random 

W(x   <    xi) 

W(x    >    xi) 

x.    =    .50    + 

l 

x.    =    1.15    - 

l 

No. 

\/w(x 

<    xi)(.2470) 

tfW(x    >    x±)  (.1755) 

.00 

.00 

(1.00) 

.50 

.01 

.01 

(.99) 

.55 

.02 

• 

.02 

• 

(.98) 

• 

.57 

• 

.58 

• 

.58 

• 

(.42) 

• 

.87 

.59 

(.59) 

.41 

.89 

.60 

• 

(.60) 

• 

.40 

• 

.89 

• 

• 

.99 

• 

(.99) 

• 

.01 

• 

1.11 

In  this  way  the  best  estimate  of  discounted  profit  per  barrel, 
corresponding  to  a  random  number  generated  to  represent  a  value 
of  W(x  <  x- ) ,  may  be  simulated.   For  use  in  the  Monte  Carlo 
Simulation,  a  Triangular  Distribution  Variable  Generator  was 
developed  and  is  included  in  Appendix  B. 

In  order  to  test  the  reliability  of  the  deviates  pro- 
vided by  the  Triangular  Distribution  Generator,  1000  values  of 
the  profit  per  barrel,  e,  were  calculated  using  the  program 
included  in  Appendix  B  and  the  values  of  a,  b,  and  c  previously 
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discussed.   The  results  of  this  test  are  listed  on  Table  7-3 
under  the  heading  of  Triangular  Distribution.   The  values  of 
profit  per  barrel  thus  simulated  were  plotted  on  Figure  7-2  as 
circles.   Also  included  on  Figure  7-2,  as  triangular  points, 
are  values  of  profit  per  barrel  calculated  manually  using  the 
technique  outlined  on  Table  7-2.   It  is  apparent  from  Figure 
7-2  that  there  is  excellent  correlation  between  the  manually- 
calculated  values  and  those  calculated  using  the  Triangular 
Distribution  Generator,  thus  confirming  the  reliability  of  the 
Generator  as  a  simulation  tool.   Dispite  this  excellent  correl- 
ation one  item  should  be  recognized.   By  the  very  nature  of  the 
random  number  generator  which  generates  values  of  numbers  used 
to  simulate  W(x  <  x-)  a  value  of  1001  for  W(x  <  x^)  can  never 
be  achieved.   Instead  the  highest  value  of  W(x  <  x^)  we  can 

achieve  using  the  generator  is  .9999 Thus,  there  is  no 

possibility  that  the  highest  value  of  profit  per  barrel  (e) 
can  ever  be  reached.   This  shortcoming  is  considered  minor 
however,  considering  the  excellent  correlation  for  virtually 
the  entire  range  of  deviates. 

7 . 6  The  Beta  Distribution  Function 

In  this  derivation  of  the  Beta  distribution  function 
Brons  [13]  notes  that  a  disadvantage  of  the  Triangular  distri- 
bution is  that  its  derivate   , ^x'    becomes  undetermined  at  x  = 

dx 

a,  b  and  c,  i.e.  the  only  points  where  some  information  about 
x  is  available.   This  means  that  at  the  outer  limits  a  and  b, 


Figure  7-2 
Cumulative  Distribution  Curve  -  Triangular  Distribution 


100 


f^i   ;;u   .u;  -+tu  -hi 

.    ;■;!  ;:::t::;;  ::::  .::; 

rrH-tHT-friT-frTTT'f  t : ;  tt , ::: 


MililsJiUln 


mpti 


LEGEND 


0  :  Computer  Generated 
A  :  Manually  Generated 


I 


50 


.60 


70 


.80        .90       1.00 
Discounted  Profit/Barrel  ($) 


1.10 


68 


w(x)  suddenly  drops  to  zero.   This  is  depicted  in  Figure  7-3 


1 

1 

1 

1 

dw(x)                              ! 
dx 

i 

Figure  7-3  -  Triangular  Distribution  Derivative 


He  also  notes  that  a  more  gentle  dropping  off  of  the  probabil- 
ity of  occurence  of  x  near  the  end  points  often  better  ex- 
presses reality  -   in  that  occurances  in  a  real  situation  would 
normally  taper  off  rather  than  end  abruptly.   This  concept  can 
be  expressed  by  the  Beta  Distribution  shown  in  Figure  7-4. 
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w(x) 


-*—  X 


Figure  7-4    Beta  Distribution  Function 

The  Beta  distribution,  like  the  Triangular  distribution,  is 
defined  by  the  three  parameters  a,  b  and  c. 
The  Beta  Distribution  is  given  by 


w 


r   _  r(a+e)xa"1(i-x)3'1 
L  J  r(cO  r(3) 


(20) 


It  can  be  seen  that  the  Beta  distribution  is  a  ratio 
of  gamma  functions  where  a  and  3  are  defined  by  the  values  of 
a,  b,  and  c.   A  derivation  of  the  relationships  between  a  and 
$  and  a,  b  and  c  (and  of  the  Beta  distribution  itself)  is  be- 
yond the  scope  of  this  thesis  and  can  be  found  in  several 
excellent  discussions  of  this  distribution  [13]  [19] . 

In  an  unpublished  report  [30],  W.G.  Lesso  developed  a 
computer  program  which  generates  deviates  having  a  Beta  dis- 
tribution.  The  program  uses  as  inputs,  the  mean  and  variance 
of  x  (as  expressed  in  a  PERT  type  analysis) : 
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(21) 


°*       36  (22) 


X 

= 

a 

+    4b 

+ 

c 

6 

2 

(c 

^2 
-a) 

From  which  the  mean  and  variance  for  the  Beta  distri- 
bution may  be  obtained  as  follows: 

x-a 

xb  "  c-a  (23) 

b  (c-a)2  (24) 

The  Lesso  developed  program,  modified  for  use  in  this  study, 
is  included  in  Appendix  C. 

For  comparision  purposes  the  Beta  Distribution  Variable 
Generator  was  used  to  calculate  1000  values  of  the  profit  per 
barrel,  e,  using  the  same  parameters  of  a,b  and  c  as  were  used 
to  calculate  the  profit  per  barrel  using  the  Triangular  Dis- 
tribution Variable  Generator.   The  values  of  profit  per  barrel 
thus  calculated  are  listed  in  Table  7-3  under  the  heading  of 
Beta  Distribution.   Cumulative  frequency  curves  for  the  Beta 
and  Triangular  distributions,  using  the  data  contained  in  Table 
7-3,  is  plotted  on  Figure  7-5.   It  can  be  seen  from  Figure  7-5 
that  very  little  difference  exists  between  the  distributions 
except  for  values  at  the  lower  end  of  the  curves.   It  should  be 
noted,  however,  that  the  mean  and  standard  deviation  for  the 
Beta  distribution  was  computed  using  the  values  as  expressed  in 
a  PERT  type  analysis  (see  Equations  21  and  22)  rather  than  in 
the  form  used  in  the  Triangular  distribution  (see  Equations  5 
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and  14) . 

Because  of  its  relative  simplicity  a  Triangular 
distribution  will  be  assumed  for  the  variable,  profit  per 
barrel.   In  the  sensitivity  studies  which  follow  however,  the 
Beta  distribution  will  be  substituted  for  the  Triangular  dis- 
tribution and  its  effect  on  the  venture  profitability  will  be 
analyzed. 
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TABLE  7-3 

COMPUTER  SIMULATED  DEVIATES  FOR  DISCOUNTED  PROFIT  PER  BARREL 

ASSUMING  TRIANGULAR  AND  BETA  DISTRIBUTIONS 

parameters 
lowest  profit  per  barrel:       $.50 
most  likely  profit  per  barrel:  $.88 
highest  profit  per  barrel:     $1.15 


TRIANGULAR 

DISTRIBUTION 

BETA  DISTRIBUTION 

Relative 

Cumulative 

Relative 

Cumulative 

rofit/bbl 

Frequency 

Frequency 

Frequency 

Frequency 

($) 

w  =   f 
1000 

f 
W   1000 

W(x  <  xi) 

W(x  <  xi) 

(%) 

W 

W 

m 

50 

0 

0 

0 

0 

51 

0 

0.0 

0 

0.0 

52 

0.2 

0.2 

0 

0.0 

53 

0 

0.2 

0 

0.0 

54 

0.5 

0.7 

0 

0.0 

55 

0.9 

1.6 

0 

0.0 

56 

0.6 

2.2 

0.1 

0.1 

57 

0.6 

2.8 

0 

0.1 

58 

0.7 

3.5 

0 

0.1 

59 

0.7 

4.2 

0.1 

0.2 

60 

0.7 

4.9 

0.3 

0.5 

61 

0.7 

5.6 

0 

0.5 

62 

0.6 

6.2 

0.6 

1.1 

63 

1.1 

7.3 

0.6 

1.7 

64 

1.3 

8.6 

0.6 

2.3 

65 

1.4 

10.0 

0.8 

3.1 

66 

1.0 

11.0 

1.0 

4.1 

67 

1.3 

12.3 

0.9 

5.0 

68 

1.2 

13.5 

1.0 

6.0 

69 

0.8 

14.3 

1.4 

7.4 

70 

1.9 

16.2 

1.2 

8.6 

71 

1.9 

18.1 

2.0 

10.6 

72 

1.4 

19.5 

2.2 

12.8 

73 

1.7 

21.2 

2.0 

14.8 

74 

1.9 

23.1 

1.6 

16.4 

75 

1.9 

25.0 

1.8 

18.2 

76 

2.0 

27.0 

2.6 

20.8 

73 


TABLE    7-3   Continued 


TRIANGULAR 

DISTRIBUTION 

BETA  DISTRIBUTION 

Relative 

Cumulative 

Relative 

Cumulative 

Profit/bbl 

Frequency 

Frequency 

Frequency 

Frequency 

($) 

f 

w  =  

1000 

£ 
w   TWO" 

W(x  <  x. ) 

W(x  <  xi) 

m 

m 

w 

m 

77 

1.8 

28.8 

2.9 

23.7 

78 

2.5 

31.3 

1.6 

25.3 

79 

2.5 

33.8 

3.4 

28.7 

80 

2.1 

35.9 

2.3 

31.0 

81 

2.3 

38.3 

2.8 

33.8 

82 

2.5 

40.7 

3.1 

36.9 

83 

3.0 

43.7 

2.3 

39.2 

84 

3.2 

46.9 

3.2 

42.4 

85 

2.6 

49.3 

3.1 

45.5 

86 

2.8 

52.3 

3.4 

48.9 

87 

3.2 

55.5 

3.7 

52.6 

88 

3.9 

59.4 

3.6 

56.2 

89 

2.6 

62.0 

3.3 

59.5 

90 

2.9 

64.9 

4.3 

63.8 

91 

2.6 

67.5 

3.7 

67.5 

92 

1.9 

69.4 

3.2 

70.7 

93 

3.2 

72.6 

2.1 

72.8 

94 

2.5 

75.1 

3.5 

76.3 

95 

2.5 

77.6 

2.1 

78.4 

96 

1.8 

79.4 

3.3 

81.7 

97 

2.5 

81.9 

2.3 

84.0 

98 

1.7 

83.6 

3.1 

87.1 

99 

2.2 

85.8 

1.6 

88.7 

100 

2.2 

88.0 

2.1 

90.8 

101 

1.5 

89.5 

1.4 

92.2 

102 

2.0 

91.5 

1.5 

93.7 

103 

1.3 

92.8 

1.3 

95.0 

104 

1.4 

94.2 

1.2 

96.2 

105 

0.8 

95.0 

0.6 

96.8 

106 

1.5 

96.5 

0.7 

97.5 

107 

0.6 

97.1 

1.0 

98.5 

108 

0.9 

98.0 

0.5 

99.0 

109 

0.4 

98.4 

0.3 

99.3 

110 

0.5 

98.9 

0.2 

99.5 

74 


TABLE  7-3  Continued 


TRIANGULAR  DISTRIBUTION 


Profit/bbl 
($) 


Relative 
Frequency 

f 
1000 

W 


w 


111 
112 
113 
114 
115 


Cumulative 
Frequency 

m 


99.1 

99.6 

100.0 

100.0 

100.0 


BETA  DISTRIBUTION 


Relative 
Frequency 


w  = 


1000 


w 


0.3 
0 

0.2 
0 
0 


Cumulative 
Frequency 

W(x  <  x±) 


99 

99 

100 

100 

100 


Figure  7-5 
Triangular  and  Beta  Distribution  Curves 
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CHAPTER  VIII 
DRY  HOLE  COSTS 

In  Chapter  V  we  discussed  the  concept  of  probability 
of  success  and  how  it  related  to  the  number  of  successful  wells 
to  be  expected  in  a  drilling  program.   To  account  for  the  num- 
ber of  wells  that  are  unsuccessful  we  have  included  dry  hole, 
or  unsuccessful  well,  costs  into  our  profitability  model. 

While  it  can  be  argued  that  dry  hole  costs  should  have 
a  continuous  frequency  distribution,  just  as  the  other  parameters 
in  the  model,  we  will  assume  for  simplicity  sake  that  the  dry 
hole  costs  are  constant  in  this  thesis. 

To  be  consistent  with  the  values  developed  in  Table 
7-1,  we  will  assume  that  the  drilling  venture  is  to  take  place 
in  the  Southern  Louisiana  on-shore  area.   From  January  to  Dec- 
ember 1970  there  were  a  total  of  1245  wells  drilled  in  this 
area  with  a  total  footage  of  12,703,309  [51].   Therefore,  the 
average  well  depth  would  be  equal  to  10,200  feet/well.   Using 
1969  average  dry  hole  costs  of  $13.23/ft  [50],  the  average  dry 
hole  cost  per  well  would  be  $135,000.   Allowing  for  inflation 
and  deeper  wells  we  will  assume  a  constant  dry  hole  cost  of 
$150,000  for  each  unsuccessful  well  in  the  program. 
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CHAPTER  IX 
COMPUTER  PROGRAM 

The  FORTRAN  IV  computer  program,  as  contained  herein, 
was  developed  to  apply  the  Monte  Carlo  simulation  technique  to 
an  exploratory  well  drilling  venture  using  the  concepts  and 
relationships  established  in  previous  chapters,  to  determine 
the  best  estimate  of  discounted  profit  and  to  accumulate  data 
for  each  iteration  so  that  a  cumulative  probability  curve  can 
be  plotted. 

The  program  is  designed  to  call  up  the  same  series  of 
random  numbers  with  each  separate  computer  run  so  that  sen- 
sitivity studies  can  be  made.   The  random  number  generator 
used  in  this  program  is  a  function  designed  for  the  CDC  6600 
computer  system  and  is  called  by  RANF(O).   The  generator  pro- 
duces pseudo  random  numbers,  distributed  uniformly  in  the  in- 
terval 0<x<l. 

The  program  accepts  as  data:   the  number  of  wells 
drilled,  probability  of  success,  dry  hole  costs,  and  the  value 
M,  which  is  the  number  of  end-points  entered  in  Line  5  plus 
one,  (Line  4);  the  median  value  of  Q,  the  value  of  Q  at  the 
median  +  1  standard  deviation,  Q84,  and  the  boundary  conditions, 
Q  max.  and  Qmin.,  (Line  6);  and  the  lowest,  most  likely  and 
highest  profit  per  barrel,  A,  B  and  C,  (Line  7).   In  addition, 
the  end  points  of  the  random  number  ranges  corresponding  to  the 
probability  of  the  number  of  successful  wells  is  entered  as  NEP, 
(Line  5) . 
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With  this  data  entered,  values  are  initialized  (Lines 
12,  13,  14,  15,  18,  20)  and  the  Standard  Deviation  and  Mean  of 
the  Lognormal  distribution  are  calculated  (Lines  16  and  17). 

For  each  iteration  specified  (Line  23)  the  first  series 
of  operations  establishes  the  number  of  successful  wells  (Lines 
24-32).   To  do  this  a  random  number  is  called,  RANF(O),  mult- 
iplied by  10,000  and  compared  with  the  end  points  corresponding 
to  the  number  of  successes  previously  entered  (Line  5).   Thus, 
the  random  number  called  establishes  the  number  of  successes, 
NSW. 

If  the  number  of  successes  is  greater  than  0,  the  pro- 
gram is  designed  to  calculate  the  size  reserve  for  each  success 
and  to  total  them  to  obtain  the  value  QSUM,  which  is  the  total 
reserve  size  discovered  (Lines  37-49) .   For  a  detailed  descrip- 
tion of  the  Lognormal  distribution  generator  see  Appendix  A. 

With  the  total  reserve  size  determined  the  profit  per 
barrel,  for  this  iteration,  E,  assuming  a  triangular  distribu- 
tion, is  calculated  (Lines  51-57).   For  a  detailed  discription 
of  the  Triangular  distribution  ganerator  see  Appendix  B. 

Having  calculated  or  inputed  all  the  parameters  the 
iteration  profit,  PROFIT,  is  computed  (Line  58)  using  the  re- 
lationship expressed  in  Chapter  III,  Equation  4.   Lines  59  and 

60  accumulate  this  iteration  profit  and  average  the  value,  PTOT, 
over  the  number  of  iterations  run,  MM.   XSQ  and  XSQTOT  in  Lines 

61  and  62  are  used  to  calculate  the  standard  deviation  and  coe- 
fficient of  variation  later  in  the  program. 
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Lines  64  through  176  are  inserted  to  store  the  values 
of  iteration  profit  in  arrays  by  PROFIT  intervals  in  order  to 
compute  the  probability  of  the  interval,  PROB,  and  the  cumu- 
lative probability,  CUMP  (Lines  189-202). 

The  variance,  VAR,  standard  deviation,  SD ,  and  coeff- 
icient of  variation,  CV,  are  calculated  for  the  particular  set 
of  model  parameters  in  Lines  178,  179  and  180  using  information 
previously  determined  in  Lines  61  and  62. 

A  listing  of  the  FORTRAN  IV  computer  program  is  con- 
tained on  the  following  pages. 


K=1,M) 
Q84,  UMAX* 


QMEDIAN,  Q84,  QMAX,  QMIN 


PROGRAM  PROFIT  ( INPUT, OUTPUT ) 

DIMENSION  N£P(25)»  IY(40>»  QLOGN(20>»  D(40>»  YY(40>»  PROB(40), 
1  CUMP(40) 
READ  100*  N»P»CU»M 
READ  110*  (NEP(K), 

READ  111,  QMEDIAN,  Q84,  QMAX,  QMIN 
READ  112.  A,  B»  C 
PRINT  155 
PRINT  113 

PRINT  114,  N»  P,  CU, 
PRINT  160,  A,  B»  C 
PROF  =  0.0 
20   PROFIT  =  0.0 
XSQ  =  0.0 
XSQTOT  =  0.0 

SD  =  ALOG(Q84)  -  ALOG (QMEDl AN) 
TMU  =  ALOG(QMEDIAN) 
MM  =  1 
PRINT  115 
DO  300  I  =  1*40 
IY(I)  =  0.0 
300   CONTINUE 

DO  99  II  =  1.1000 

CALL  RANF(O) 

Ml  =  RANF(O) 

w  =  WlMOOOO. 

J  =  w 

DO  10  I  =  l.M 

IF  (J.Le.NeP(D)  GO  TO  11 

10  CONTINUE 

11  NSW  =  1-1 
IF  (NSW. FQ. 0.0)  32*40 

32   QSUM  =0.0 
E  =  0.0 
GO  TO  56 

40   CONTINUE 

QSUM  =  0.0 

DO  45  I  =  l.NSW 

RA  =  RANF(O) 

RB  =  RAnF(O) 

V  =  <-2.0»ALOG(RA))«*.5*COS(2.0«3.1416*RB) 

QNORM  =  V»SD  ♦  TMU 

QLOGN(I)  =  EXP (QNORM) 

IF  (QLOGN(I)-QMIN)  46*47,48 

46  QLOGN(I)  =  QMIN 
GO  TO  47 

48  IF  (QLOGN(I)-QMAX)  47,47,49 

49  QLOGN(I)  =  QMAX 

47  QSUM  =  QSUM  ♦  QLOGN(I) 
45   CONTINUE 

CALL  RANF(O) 

WXLTXI  =  RANF(O) 

WXGTXI  =  1.0  -  WXLTXI 

XI  =  A  ♦  SQRT(WXLTXIMB-A)«(C-A)) 

IF  (XI.LE.B)  55,50 

50  XI  =  C  -  SQRT(WXGTXI*(C-B)«(C-A)) 

55  E  =  XI/100. 

56  PROFIT  =  QSUM«E  -  (N-NSW)«CU 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 
26 
27 

28 
29 

30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 


bl 


62 


63 


64 


65 


66 


67 


PnOF  =  PROF  I T  ♦  PROF 

PTOT  =  P&OF/MM 

XSQ  =  (PROFIT)**2 

XSQTOT  =  XSQTOT  ♦  XSQ 

PRINT  116*  MM»  NSW» 

QSUM* 

F,  PROFIT. 

mm  =  mm  ♦  i 

IF(PROFIT.LT.-IO.O) 

GO 

TO 

61 

IFtPROFIT.LT.  0.00) 

GO 

TO 

62 

IF (PROFIT. LT.  10.0) 

GO 

TO 

63 

IF(PR0FIT,LT.  20.0) 

GO 

TO 

64 

IF  (PROFIT. LT.  30.0) 

GO 

TO 

65 

IF  (PROFIT. LT.  40.0) 

GO 

TO 

66 

IF (PROFIT. LT.  50.0) 

GO 

TO 

67 

IF (PROFIT. LT.  60.0) 

GO 

TO 

68 

IF(PROF 1T.LT.  70.0) 

GO 

TO 

69 

IFtPROFIT.LT.  80.0) 

GO 

TO 

70 

IF(PROFIT.LT.  90.0) 

GO 

TO 

71 

IF (PROFIT. LT. 100.0) 

GC 

TO 

72 

IF  (PROFIT.LT. 110.0) 

GC 

TO 

73 

IFtPROFlT.Lr. 120.0) 

GC 

TO 

74 

IF(PROFiT.LT. 130.0) 

GO 

TO 

75 

IF(PROFIT.LT. 140.0) 

GO 

TO 

76 

IF (PROFIT.LT. 150.0) 

GO 

TO 

77 

IF(PROFIT.LT. 160.0) 

GO 

TO 

78 

IF<PROf IT.LT. 170.0) 

GO 

TO 

79 

IF (PROFIT.LT. 180.0) 

GO 

TO 

80 

IF (PROFIT.LT. 190.0) 

GO 

TO 

81 

IF (PROFIT.LT. 200.0) 

GO 

TO 

82 

IF (PROFIT.LT. 210.0) 

GO 

10 

83 

IF (PROFIT.LT. 220.0) 

GO 

TO 

84 

IF (PROFIT.LT. 230.0) 

GO 

TO 

85 

IF(PROFIT.LT. 240.0) 

GO 

TO 

86 

IF(PROFIT.LT. 250.0) 

GO 

TO 

87 

IF (PROFIT.LT. 260.0) 

GO 

TO 

88 

IF (PROF IT.LT. 270.0) 

GO 

TO 

89 

IF (PROFIT.LT. 280.0) 

GO 

TO 

90 

IF (PROF IT.LT. 290.0) 

GO 

TO 

91 

IF  (PROFIT.LT. 300.0) 

GO 

TO 

9? 

IF<PR0FIT.LT. 350.0) 

GO 

TO 

93 

IF(PR0FIT.LT. 400.0) 

GO 

TO 

94 

IF(PROF IT.LT. 450.0) 

GO 

TO 

95 

IF (PROFIT.LT. 500.0) 

GO 

TO 

96 

IF(PROFIT.LT. 550.0) 

GO 

TO 

97 

GO  TO  99 

IY  (1)  =  I Y (  1  )  ♦  1 

GO  TO  99 

IY(2)  =  IY(2)  ♦  1 

GO  TO  99 

IY (3)  =  IY(3)  ♦  1 

GO  TO  99 

1Y(4)  =  IY<4)  ♦  1 

GC  TO  99 

IY(5)  =  1Y(5)  ♦  1 

GO  TO  99 

IY(6)  =  IY(6)  ♦  1 

GO  TO  99 

IY(7)  =  IY(7)  ♦  1 

GC  TO  99 

PTOT 


59 

60 
61 
62 
6  3 
6w 
bb 
66 
b7 
68 
69 
70 
71 
7? 
73 
74 
7b 
7b 
77 
78 
79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

116 


68  I Y ( 8  >  =  IY(B)  ♦  1 

(30  TO  99 

69  IT(9)  =  IY(9)  ♦  1 
GO  TO  99 

70  IYUO)  =  IY  (10 
GO  TO  99 

71  ir (11 )  =  I y  <  1 1 

GO  10  99 

72  IY ( 12)  =  I Y ( 12: 

GO  1 0  99 

73  IY  (13)  =  I Y  < 1 31 

GO  TO  99 
Ik        I  Y  (  1**)   =  I  Y  (  14: 

GO  TO  99 
7^   I  Y  (  15)  =  I  Y  (15! 

GO  10  99 

76  I Y ( lb)  =  I Y (16: 
GO  TO  99 

77  1YU7)  =  I  Y  (  1  7 
GO  10  99 

fti        1 Y (18)  =  I Y (18: 
GO  10  99 

79  I Y ( 19)  =  I T (19: 
GO  10  99 

80  IY<20)  =  IY(20: 
GO  TO  99 

HI   IY (21  )  =  IY (21 

-C  10  -*9 
H2   1Y(22>  =  IY(22 

■jO  10  99 
H3   IY(23)  =  IY(23: 

ijG  1U  99 

84  Iy(24)  =  IY(24: 
GO  10  99 

85  IY(25>  =  IY(25: 
GC  TO  99 

se   I Y (26)  =  IY (26: 
uO  TO  99 

87  IY(27)  =  IY(27: 
Go  10  99 

88  IY(28)  =  IY(26: 
GC    TO  99 

89  IY(,>9)  =  IY(29: 
GO  TO  99 

90  IY(30)  =  IY(30: 
GO  TO  99 

91  IY (31 >  =  IY (31 
GO  TO  99 

92  IY(32)  =  IY(32: 

GO  TO  99 

^3       IY<33)  =  IY(33: 

GO  TO  99 
94   IY<34>  =  IY(34: 

GO  TO  99 
s5   IY(35)  =  IY(35: 

GO  TO  99 
96   IY<36)  =  IY(36: 

GO  TO  99 


117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 

12* 
129 
13o 
131 
132 
133 
134 
135 
136 
137 
138 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
151 
152 

153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
1*4 
165 
166 
167 
168 
169 
170 
171 
172 
173 
174 
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PRINT  140* 

MMfPTOT 

PRINT  14 1, 

VARt  SDt  CV 

PRINT  144 

PRINT  155 

PRINT  113 

PRINT  114» 

N»  Pf  CU.  QMEDIAN 

PRINT  160» 

A»  B»  C 

PRIM  145 

SUMM  =  0.0 

97   IY(37)  =  IY(37)  ♦  1 

99   CONTINUE  176 

60   MM  =  MM  —  1  177 

VAR  =  (XSQTOT  -  (MM* <PT0T**2) ) ) / (MM  -  1)  178 

SO  =  SURT(VAR)  179 

CV  =  SD/PTUT  180 

181 
182 
183 
184 
185 

l»  Q84»  UMAX,  QMIN  186 

187 
188 
189 

0(1)  =  -10.0  190 

00  200  IM  =1.37  191 

YY(IM)  =  IY(IM)  192 

PROMIM)  =  YY(IM)/1000.  193 

SOMM  =  SUMM  ♦  PROm(IM)  194 

CUMP(IM)  =  SUMM  195 

KJ  =  IM  ♦  1  196 

IF(IM-32)  210.  220.  220  197 

210   O(KJ)  =  Q( IM)  ♦  10.0  198 

GO  TO  230  199 

220   D(KJ)  =  D(IM)  ♦  50.0  200 

230   CONTINUE  201 

PHINT  150.  D(IM),  PROm(IM),  CUMP(IM)  202 

200   CONTINUE  203 

100   FORMAT  (10X,I5»2F10.3.I5)  204 

110  FORMAT  (10X.17I4)  205 

111  FORMAT  (10X,4F10.3)  206 

112  FORMAT  (10X,3F10.2)  207 

113  FORMAT  <///,6X,*N0.  OF  WELLS*»3X ,*PROB  0F«»2X»»DRY  HOLE**/ t9X,»DR I  208 
1LLED**5X,*SUCCESS*»4X,*C0ST**4X,»QMEDIAN*»3X,*Q84*»5X»*QMAX*,5X,  209 
2*0MlN«)  210 

114  FORMAT  <10XtI5»8X,F4.3.4XfF4.3  OX»F7.2*2X»F7.2*2X»F7.2»2X»F7 .2//)  211 

115  F0RMAT(5X,*ITERATlON*t3X,»NUMBER  0F*,4X,*FIELD*,2X t *PROF IT  PER*t2X  212 
l.*ITE«ATlON«,3Xt*BEST  ESTIMATE* »/l 7X. "SUCCESSES* *4X,«S1ZES*.5X »»BA  213 
2RREL»»6X,«PR0FIT»»6X»*0F  PROFIT*,//)  214 

116  FORMAT  (   5XfI4tllX»I3   ♦6X»F6.2 ♦6X.F4.2* 7X»F6.2»8X ,F6.2>  215 

140  FORMAT  (///♦SX^THE  BEST  ESTIMATE  OF  PROFIT  BY  THE  MONTE  CARLO  SIM  216 
1ULATION  METHOD*»/7X,*AT  *»I4<*  ITERATIONS  IS  *,F7.2,*  MILLION  DOLL  217 
1AkS«>  218 

141  FORMAT  (///♦5X,»THE  VARIANCE  IS  *»F12.3t/»5X, *THE  STANDARD  DEVIATI  219 
ION  IS  *,F12.3»/»5X,*THE  COEFFICIENT  OF  VARIATION  IS  *,F6.4>  220 

14^   FORMAT  (lHl)  221 

145   FORMAT     (15Xt*UPPER  LIMIT  OF*»7X,*PROBABIL ITY*»5X,«CUMMULAT IVE* ♦  222 

1/»15X,*PR0FIT  INTERVAL*. 5X»*0F  INTERvAL*»5X.»PR0BABILITY*,///>  223 

150   FORMAT  (,?0X,F4.0»15X.F5.4.11X.F5.4>  224 

155   FORMAT  < lHl ♦ 15X,»PR0F  IT  DISTRIBUTION  OF  AN  EXPLORATORY  DRILLING  PR  225 

10GRAM*)  226 

160   FORMAT  (/15X,*THE  TRIANGULAR  DISTRIBUTION  FOR  DETERMINATION  OF  THE  227 

1*,/15X,*PR0FIT  PER  BARREL  HAS  THE  FOLLOWING  EXTREMES  IN  CENTS  -  *.  228 

2/»l5Xt*A  =  *tF6.2»7X»*B  =  *,F6.2,7X»*C  =  *»F6.2»///>  229 

EinD  230 


CHAPTER  X 
BEST  ESTIMATE  OF  DISCOUNTED  EXPECTED  PROFIT  USING  BASE  DATA 

In  previous  chapters  the  profitability  model  was  dev- 
eloped and  the  individual  parameters  and  their  frequency  dis- 
tributions were  discussed.   In  addition,  certain  Base  Data  was 
established  for  each  parameter  so  that  the  best  estimate  of 
discounted  expected  profit,  E(P),  could  be  calculated.   This 
Base  Data  is  summarized  in  Table  10-1. 

TABLE  10-1 
BASE  DATA 


Parameter 

Distribution 

Value 

Number  of  Wells  Drilled  -  n 

- 

20 

Probability  of  Success  -  p 

Binomial 

.15 

Reservoir  Size 

Lognormal 

Q  at  51  probability 

5 

X 

106  bbl 

Q  at  951  probability 

79 

X 

106  bbl 

Q  maximum 

140 

X 

106  bbl 

Q  minimum 

1 

X 

106  bbl 

Discounted  Profit/Barrel 

Triangular 

Lowest  -  e0 

$ 

.50/bbl 

Most  Likely  -  e^ 

$ 

.88/bbl 

Highest  -  ec 

$] 

..15/bbl 

Dry  Hole  Cost  -  Cu 

None 

$ 

150,000 

Number  of  Iterations 

1,000 
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On  the  basis  of  this  Base  Data  and  using  the  Fortran 
IV  program  included  in  Chapter  IX  a  value  of  Best  Estimate  of 
Discounted  Expected  Profit,  E(P),  of  $68.02  million  was  obtained 
See  Appendix  D  for  a  printout  of  the  number  of  successes, 
reserve  sizes,  profit  per  barrel,  iteration  profit  and  average 
best  estimate  of  discounted  profit  for  each  of  the  1000  iter- 
ations using  the  Base  Data.   In  addition,  Appendix  D  contains 
cumulative  probability  data  for  each  specified  interval  of  E(P). 
Using  this  cumulative  probability  information  the  histogram 
included  as  Figure  10-1  was  prepared  from  which  it  can  be  seen 
that: 

1.  The  Median  Value  (at  50%  probability)  of  Discounted 
Profit  is  $56.03  million. 

2.  With  90%  confidence  we  could  expect  a  discounted 
profit  of  between  $1.51  million  and  $170.0  million. 

3.  The  probability  of  a  loss  is  4.2%,  or  stating  it 
another  way,  there  is  a  95.8%  chance  of  a  profit. 

Further,  from  the  printout  in  Appendix  D  we  see  that: 

4.  The  standard  deviation  is  $53.66  million. 

5.  The  coefficient  of  variation  is  .7889. 

The  significance  of  the  standard  deviation  and  coefficient  of 
variation  as  measures  of  risk  were  discussed  in  Chapter  IV. 


Figure  10-1 
Cumulative  Probability  -  Base  Data 
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CHAPTER  XI 
SENSITIVITY  ANALYSIS 

Freedom  from  uncertainty  is  a  luxury  rarely  enjoyed  by 
the  contemporary  decision  maker.   In  a  society  characterized  by 
change,  uncertainty  is  an  established  and  accepted  fact  of  life. 
This  is  especially  true  in  the  oil  business.   Notwithstanding 
the  difficulties,  the  decision  maker  is  expected  to  attain 
satisfactory  results  as  measured  by  the  goals  of  his  organiza- 
tion. 

In  the  face  of  uncertainty,  the  most  recurring  quest- 
ions he  faces  are  of  the  form,  "What  if...? 

1.  What  if  labor  costs  rise  and  dry  holes  cost  twice 
as  much?   Three  times  as  much? 

2.  What  if  the  posted  price  of  oil  decreases  by  30%? 

3.  What  if  we  drill  30  exploratory  wells  rather  than 
20?   How  about  10  wells? 

The  "what  if"  question  may  be  viewed  as  an  introduction  to  ths 
sensitivity  analysis.   In  its  mathematical  sense,  sensitivity 
analysis  is  a  study  to  determine  how  possible  changes  or  errors 
in  parameter  values  affect  model  outputs.   In  its  applied  or- 
ganizational setting,  sensitivity  analysis  may  be  broadly  de- 
fined as  a  study  to  determine  the  responsiveness  of  the  con- 
clusions of  an  analysis  to  changes  or  errors  in  parameter 
values  used  in  the  analysis.   In  other  words,  would  our  decision 
be  any  different  if  a  parameter  value,  such  as  the  probability 
of  success  in  out  profitability  model,  were  to  change.   The 

87  . 


88 


prime  function  of  sensitivity  analysis  can  be  said  to  be  to 

provide  a  better  understanding  of  risk.   Sensitivity  analysis 

can  provide  guidelines  for  allocating  organizational  resources 

to  data  collection  and  data  refinement  activities  [40]  .   The 

benefits  of  subjecting  models  to  sensitivity  analysis  are 

summarized  by  Arnoff  and  Netzorg  [2] : 

The  use  of  an  operations  research  model  is  especially 
important  and  advantageous  in  that  it  permits  experi- 
mentation "on  paper",  without  manipulation  of  the  act- 
ual system.   In  using  the  model,  one  can  assess  the 
sensitivity  (response)  of  the  system  to  a  wide  variety 
of  conditions  -  without  requiring  either  the  time, 
expense,  or  risks  associated  with  experimenting  with 
the  system  itself  (if  such  experimentation  would,  in 
fact,  be  possible  and  meaningful).   Hidden  relation- 
ships can  be  brought  to  light  and  brought  to  bear  upon 
decisions  and  control  of  activity. 

In  this  thesis  a  profitability  model  for  an  exploratory 
well  drilling  venture  was  developed  and  tested  with  Base  para- 
meters.  A  decision  problem  is  now  faced.   If  the  model  was 
developed  in  order  to  provide  results  which  would  serve  as  a 
basis  for  recommending  action  to  management,  it  is  possible, 
at  this  point,  to  use  the  model  results  as  is,  or  to  select 
other  alternatives:   to  collect  additional  data;  refine  initial 
parameters:  or  even  alter  the  model  itself.   On  the  significance 
of  this  decision,  Willian  T.  Morris  [35]  comments: 

The  nature  of  management  science  is  such  that  sooner  or 
later  it  has  to  recommend  some  action  other  than  making 
more  observations.   Usually  this  comes  sooner  rather 
than  later.   Eventually  it  will  suggest  that  action  be 
taken  consistent  with  either  acceptence  or  rejection 
of  the  hypothesis.   It  will  recommend  management  "act 
as  if"  the  hypothesis  were  true  or  "as  if"  it  were 
false.   Indeed,  the  fundamental  question  in  management 
science  may  well  be  when  to  stop  collecting  data  and 
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developing  the  model,  and  when  to  produce  a  recommen- 
dation for  action. 

In  the  past  the  emphasis  has  been  on  whether  information  re- 
quired for  decision  making  could  be  developed.   With  our  grow- 
ing capacity  to  generate  information  or  data,  the  emphasis  now 
appears  to  be  shifting  to  the  question  as  to  what  extent  this 
capacity  can  be  utilized  beneficially.   It  may  be  recognized 
that  while  the  collection  and  refinement  of  data  are  costly, 
it  may  serve  to  improve  recommentations  for  action.   In  most 
cases,  however,  little  basis  exists  for  projecting  whether  the 
additional  information  is  worth  the  cost.   Statistical  decision 
theory,  and  a  model  such  as  the  one  developed  in  this  thesis, 
will  assist  in  the  decision  as  to  whether  additional  informa- 
tion is  "worth  the  cost".   The  basic  requirements  for  the  stat- 
istical decision  theory  model  are  summarized  by  Morris  [36]  as 
follows : 

[The  Model]  must  be  able  to  enumerate  the  possible 
outcomes  of  future  data  collection  efforts  and,  further, 
to  compute  the  probabilities  of  these  outcomes.   In 
addition,  it  must  be  possible  to  indicate  just  how  the 
information  willquantitatively  change  the  decision 
maker's  view  of  his  choice. 

In  an  excellent  example,  which  illustrates  this  concept,  C. 
Jackson  Grayson  [21] ,  demonstrates  the  application  of  decision 
theory  to  drilling  decisions  by  oil  and  gas  operators  and  also 
to  the  decision  problem  of  whether  to  purchase  additional  in- 
formation before  making  a  final  drilling  decision. 

Clearly  then,  a  sensitiviey  analysis  which  discloses 
decision  insensitivity  to  estimated  parameter  variations 
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precludes  the  necessity  of  introducing  more  costly,  time 
consuming  data,  collection  data  refinement  and  statistical 
decision  models  into  the  analysis. 

It  is  in  this  vein  that  the  sensitiviey  studies  are 
conducted  in  this  thesis.   The  profitability  model  provides 
the  necessary  vehicle  for  the  analysis.   From  the  results  it 
is  possible  to  determine  the  outcomes,  as  well  as  probabilities 
of  these  outcomes,  associated  with  changes  in  parameter  values 
thus  providing  some  answers  to  the  "What  if?"  questions  posed 
at  the  beginning  of  this  chapter. 


CHAPTER  XII 
SENSITIVITY  STUDIES 

12.1  Outline  of  Studies 

The  sensitivity  studies  were  conducted  using  the 
profitability  model  developed  in  Chapter  III.   In  these  studies 
each  parameter  in  the  model  was  varied  from  the  base  values  as 
summarized  in  Table  10-1.   Each  time  a  parameter  was  varied  all 
remaining  parameters  were  set  to  their  base  values  and  the  FOR- 
TRAN IV  program,  developed  in  Chapter  IX,  was  executed  1000 
times.   In  this  manner  a  printout  similar  to  the  one  for  the 
Base  data,  included  in  Appendix  D,  was  obtained  for  each  para- 
meter change  or  Case. 

Thus,  for  each  Case,  at  1000  iterations,  we  have  the 
best  estimate  of  discounted  profit,  the  standard  deviation,  the 
coefficient  of  variation,  the  iteration  profit,  cumulative 
probability  data,  etc.   Tables  12-1  through  12-9,  grouped  by 
parameter  being  varied,  contain  a  summary  of  some  of  this  data 
(at  1000  iterations)  as  well  as  extrapolated  information  such 
as  the  median  discounted  profit,  the  maximum  and  minimum  dis- 
counted profits  within  90%  confidence  limits  and  the  chance 
of  making  a  profit  greater  than  zero.   The  values  using  the 
Base  data  have  been  included  on  each  Table  to  assist  in  eval- 
uation.  It  should  be  noted  that  each  time  the  model  was  run 
the  same  series  of  random  numbers  were  called  up  so  that  com- 
parisons between  Cases  could  be  made.   Figures  12-1  through 
12-8  contain  histograms  for  several  Cases,  plotting  the  cum- 
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ulative  probability  for  each  specified  interval  of  discounted 
expected  profit.   Figures  12-1  through  12-8  are  similar  in 
presentation  to  Figure  10-1  for  the  Base  data.   On  each  hist- 
ogram the  90%  Confidence  Limits  have  been  represented  by  dashed 
lines  at  5%    and  95%  cumulative  probabilities. 

To  assist  in  the  comparison  of  each  Case  to  the  Base, 
Tables  12-1A  through  12-9A  have  been  prepared  and  are  inserted 
behind  the  Table  it  serves  to  amplify.   These  Tables  contain 
the  percentage  differences  between  the  Case  parameter  set 
values  and  the  Base  parameter  set  values.   The  percentages 
were  calculated  using  the  following  relationship: 

A  fin  n  =  (Case  data  value)  -  (Base  data  value! 
1    0)  (Base  data  value) ~x   10°  >   U) 

In  addition  to  the  Cases  noted  above,  the  Base  data 
was  run  for  10,000  iterations  rather  than  1000  iterations. 
Additionally,  the  Beta  Distribution  Variable  Generator  was 
substituted  for  the  triangular  distribution  of  profit/barrel 
in  the  FORTRAN  IV  Program  and  the  Base  data  was  run  for  1000 
iterations.   The  results  of  both  these  variations  are  summar- 
ized on  Table  12-10. 


12.2   General  Observations 

One  of  the  prime  purposes  of  the  sensitivity  studies 
is  to  assess  the  change  in  observed  values  resulting  from  a 
change  in  parameter  value.   The  majority  of  the  effects  of 
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Figure  12-1 
Cumulative  Probability  -  'n'  Parameter  Set 
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Figure  12-2 
Cumulative  Probability  - 'p ' Parameter  Set 
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Figure  12-3 
Cumulative  Probability  -  Q  at  5%  Parameter  Set 
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Figure  12-4 
Cumulative  Probability  -  Q  at  95%  Parameter  Set 
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Figure  12-5 
Cumulative  Probability  -  Cu  Parameter  Set 
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Cumulative  Probability  -  e(a)  Parameter  Set 
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118 


T" T'TT " " "  j  "mmmJ9 ' 

s  — _._ 1.   ___± ::::::::::::::: -y-  - 

°  —  v- |._-i___i t           X — ?9  i ! 

s i_u. case  53    yA::::::;^ 

:^:::z:"=||  Ji=:::::::z:EEEEEEEEEEE::     $  •  62/bbl    ::-^--::::-::5--:± 

-t----f--        r'"tT" >^     1 -*-— F-*-*- --  j 

'''                  ^       ' 

/'                    _** 

> '                i>^ 

sli|44j^                                          |  1  j  1 1 1  ■  1 1 1 ' '  i '  | '  1 ' '  1 1  1 1  j 

1— i H ' —J i -»- ■     >         1   i    i                                          1                                 -l-                — —    -1 .       i 

■  -J--*- hi--1-1 — *-   H  '  i H b*- 44-— -*-r-      Mi U-t- U ,lii     ,  i  i      i  ! 

~m~ !  -"r-4-  |  \y 7rr ' H 1 TTTT;       ~m~  1  !  !  !  "Tr "  !  1  ;  l~  i  !  !  ! 

i                              i       i     i     /\        S                                                                       i 

_|_          i  1  Ml/ r  \     /            1  i                                         'II                  ill 

i                                        Ml           III/'      1   |    y'                    i| 

o  ■                  i  \  \\\y     ]/                     i                   '           -                    i 

™  Tr        i f     \\jr      [if 

1       ■                  '  X  ■    l/i                         1 

1                              !         i/l  i    jf?                                                                                             ...........                                j- 

/  i/i                !ti       M               ^  i^  i      ,    '                           t 

5  -1--U.-U-I J/|li>flillllllllllll  ^-^Li. L  +  --u-a-il L_.L-._4- 

-  ol|==#=Pll   S?5?^  :iSEEEglEE,EEEEzEEb_EEEEEEcEEE_E*E*± 

*     SE^Ep^^    $1.14/bbl:^#i#==^^±=fe=^Kf 

AR     o t-t--1 h'^-+/+--h — M — H 1 1  i  1  i  -1 — H--\ 1 1 ! >+-\ L 

g     MM \J  1  \f\  .      Mil      III _U-  -III    — 1—  —  -J — 1 — 1 L. -i U4_    M  |  , 

*                     1                /     A  II                                                                                                                                                                        1 

n                     1                /    /                                                                                                                                                                                        ' 

w              1    :     Jfl  r                          1                                                                                                , 

"*■        ■*        M             /  l/           1      i                  II                              '                     1      i                                '                                                          ' 

"■                fir     ft                                                      m(          j    i                                  [  1 

/    Jr 

J     /I 

g  -XL-tj-l H+ ' M 1 1 — 1 — H ^ — ^"^ 1 

_    _ 1j— i—l—      _ 1— _ i__j ;_^_  _      __      — i— — 1 ■_)___. u 

Jy jl''  hi i  ' hi     1       '      11'    11     hi  11  i 111    111    Ml  Nil 
° -Jl              '  i                        \    —  ll. I__l___E___ 

(/                                        i  1                                                                                                                                                    

I(                                                                                                                                                                                       1                   1 

■  __U  =  l-±-=--it=U^ 

i H     '     '  H mi                           I       '            r   i        M 

i_                                 i 

«  _J 1  i  1  i  -L.-U J- J 1 — L-I L 1 L-U-  1  |  U L--L 

<*      •  ■  1     1  ■         1  H  I       I  '  :                '  1  1       [  •  II     II  |       1         lift  1 1             ;i     ;                        T 

___p:_t=__z_±4::__t1:±p-_=|=q= 1 •--- f-H ' ' 

EE::      ~~    % _t__ ::___::__  ::_  _::___::_::::_—::    ::::_:::::± 

^^_^.j^___                                                     h    __ 

jl IT _t___+ 

■> 

S                                                   TL                   " 

__L_±±-L 

0                 kO               80             120             160             200             2*0             280      ' 

X 


PROFIT  (In  MllUon  Dollars) 


Figure  12-8 
Cumulative  Probability  -  e(c)  Parameter  Set 
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these  parameter  changes  can  be  observed  on  Tables  12-1  through 
12-9  and  12-1A  through  12-9A  and  require  no  further  datailed 
discussion.   The  analysis  of  a  parameter  change  to  a  change  in 
Best  Estimate  of  Discounted  Profit  and  to  the  risk  as  measured 
by  the  changes  in  Coefficient  of  Variation  will  be  discussed 
in  the  next  Section. 

At  this  point  however,  some  general  observations 
regarding  the  sensitivity  of  the  model  to  parameter  changes 
can  be  made.   First,  an  inspection  or  the  Tables  will  reveal 
that  the  one  single  observed  value  which  fluctuated  the  great- 
est on  a  percentage  basis  was  the  minimum  value  of  discounted 
profit  which  could  be  expected  with  90%  confidence.   The  reason 
for  this  high  percentage  fluctuation  was  the  relatively  small 
numbers  involved.   Thus,  a  variation  of  $7.87  million  on  a 
$1.51  million  Base  resulted  in  a  percentage  change  of  approx- 
imately 980%. 

The  second  observetion  made  from  the  Tables  is  that, 
except  in  isolated  cases,  the  chance  of  a  profit  greater  than 
zero  was  fairly  constant  at  approximately  95  to  961.   Stating 
it  another  way,  there  appears  to  be  only  a  4  to  51  chance  of  a 
loss  in  a  majority  of  the  Cases.   When  the  probability  of 
success  was  reduced  to  1  in  20  there  was  still  a  621  chance  of 
making  a  profit  and,  even  when  the  number  of  wells  drilled  was 
reduced  to  5  rather  than  20,  the  probability  of  a  profit  was 
approximately  561  (the  lowest  value  in  any  Case) . 

An  inspection  of  Table  12-10  reveals  that  the  substitu- 
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tion  of  the  Beta  distribution  for  discounted  profit  per  barrel 
rather  than  the  Triangular  distribution  has  very  little  effect 
on  the  Best  Estimate  of  Discounted  Profit,  E(P).   The  risk  as 
measured  by  the  Standard  Deviation  and  the  Coefficient  of 
Variation  is  increased  however,  by  about  4.41.   The  minimum 
value  of  discounted  profit  with  901  Confidence  Limits  has  been 
decreased  by  4.341  and  the  maximum  value  has  been  increased  by 
1.78%.   Additionally,  the  range  of  values  with  90%  Confidence 
Limits  has  been  increased  with  the  Beta  distribution  by  2.24%. 
With  regard  to  the  Base  data  run  at  10,000  iterations 
rather  than  1000  iterations,  it  is  seen  from  Table  12-10  that 
the  value  of  E(P)  only  varies  by  1.47%  from  the  1000  iteration 
value  while  the  relative  risk  as  measured  by  the  Coefficient 
of  Variation  increases  by  1.5%  primarily  due  to  the  decrease 
in  the  absolute  value  of  E(P). 

12.3  Effects  of  Parameter  Changes  on  the  Best  Estimate  of 
Discounted  Profit  and  Coefficient  of  Variation 

The  first  question  explored  in  this  Section  is  the 
effect  of  a  change  in  parameter  value  on  the  best  estimate  of 
discounted  profit,  E(P).   This  effect  shall  be  assessed  by 
calculating  a  Significance  Ratio  for  each  Case.   The  Signifi- 
cance Ratio  is  defined  as  the  change  in  observed  value  (i.e. 
E(P))  to  a  change  in  the  parameter  value.   In  a  conceptual 
sense  the  Significance  Ratio  measures  the  significance  of  a 
change  in  a  parameter  of  the  model  on  a  value  calculated  by 
the  model.   We  assume,  of  course,  that  all  other  parameters 
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(except  the  one  being  studied)  remain  the  same.   Mathematically, 
this  can  be  expressed  as: 

A  Observed  Value 


SR  = 


5  Parameter 


(1) 


With  this  Significance  Ratio  the  relative  importance  of  a 
change  in  various  model  parameters  on  the  value  being  observed 
can  be  measured.   For  the  profitability  model  it  is  possible, 
for  example,  to  measure  the  relative  importance  of  a  change 
in  dry  hole  costs,  lowest  profit  per  barrel,  etc.  on  the  best 
estimate  of  discounted  profit  for  the  total  venture.   For  the 
analysis  of  the  best  estimate  of  discounted  profit,  Equation 
1  can  be  rewritten: 

AE(P)  


SR  = 


A  Parameter 


(2) 


The  values  to  be  inserted  in  Equation  2  are  obtained  from 
Tables  12-1A  through  12-9A.   As  an  example,  from  Table  12-1A 
the  following  Significance  Ratio  is  obtained: 

For  a  change  in  parameter  of  -75.00%  (Case  1)  from  the 
Base  we  have  a  commenserate  reduction  of  E(P)  of  -76.77%. 
Inserting  these  values  in  Equation  2  we  have: 


SR  = 


-76.77 
-75.00 


SR  =  1.02 
For  Case  2  we  have  the  following  relationship 

-50.44 
SR  "    -50.00 

SR  =  1.0. 
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In  a  similar  fashion  the  Significance  Ratios  can  be  calculated 
for  each  of  the  nine  Cases  involving  a  variation  of  the  number 
of  wells  drilled  as  shown  on  Table  12-1  and  12-1A.   If  this  is 
done  Significance  Ratios  ranging  from  .81  to  1.18  with  a  mean 
value  of  0.99  are  obtained. 

Table  12-11  contains  a  summary  of  the  Ranges  and  Mean 
Significance  Ratios  for  each  of  the  parameter  sets  contained  on 
Tables  12-1  through  12-9   and  Tables  12-1A  through  12-9A.   In 
addition,  Table  12-11  contains  a  value  of  the  Normalized 
Significance  Ratio  which  reduces  each  Significance  Ratio  to  a 
value  on  a  1-100  scale  in  proportion  to  its  actual  Significance 
Ratio.   The  Normalized  Significance  Ratio  can  be  utilized  to 
facilitate  comparison  of  parameter  sets 

TABLE  12-11 
SIGNIFICANCE  RATIOS  -  E(P) 

Significance  Ratios 
Parameter  Set         Range       Mean       Normalized  Mean 

96 
100 
29 
70 
07 
0 
04 
20 
35 
46 


n 

.81 

-    1 

.18 

.99 

P 

.96 

-    1 

.10 

1.03 

Q   at    5% 

.19 

-    • 

40 

.30 

Q   at    951 

.56 

-    . 

82 

.72 

Q(max) 

.02 

"    • 

12 

.07 

Q(min) 

0 

0 

Cu 

0 

.04 

e(a) 

0 

.21 

e(b) 

0 

.36 

e(c) 

0 

.47 
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Several  observations  can  be  made  from  Table  12-11. 
First,  it  can  be  seen  that  a  change  in  the  number  of  wells 
drilled  and  probability  of  success  are  both  almost  equally 
significant  as  far  as  their  effects  on  the  value  of  E(P).   In 
fact,  there  is  virtually  a  1  to  1  relationship  between  a  change 
in  the  number  of  wells  drilled  and  probability  of  success  and  a 
change  in  the  best  estimate  of  discounted  profit.   In  other 
words,  if  the  number  of  wells  drilled  is  increased  by  10%, 
almost  a  101  increase  in  discounted  profit  can  be  expected. 
Likewise,  if  the  probability  of  success  is  decreased  by  201 
about  a  20%  decrease  in  the  discounted  profit  can  be  expected. 

Also  seen  from  Table  12-11  is  that  the  next  most 
significant  parameter  in  the  model  is  the  estimated  value  of 
reserve  size  with  95%  probability,  Q  at  95%.   In  fact,  this 
parameter  has  a  normalized  significance  of  70  which  is  approx- 
imately 2.5  times  as  critical  as  the  lower  estimate  of  the 
reserve  size,  Q  at  5%  probability,  on  the  discounted  profit. 
In  a  similar  fashion  it  can  be  observed  that  variations  in  the 
values  of  the  upper  and  lower  boundary  limits  for  reserve  sizes 
have  little  or  no  significance  on  the  discounted  profit. 

Although  one  might  intuitively  feel  that  changes  in  dry 
hole  costs  would  have  a  material  effect  on  the  discounted 
profit  it  is  seen  that  it  only  has  a  normalized  significance 
of  4  which,  in  essence,  is  relatively  insignificant  in  compar- 
ison to  the  other  parameters. 

In  comparing  the  three  parameters  associated  with  the 
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discounted  profit  per  barrel  it  is  noted  that,  of  the  three, 
changes  in  the  highest  profit  per  barrel,  e(c),  are  2.3  times 
as  significant  as  changes  in  the  lowest  value,  e(a),  and  1.3 
times  as  significant  as  changes  in  the  most  likely  value,  e(b) 
Also  observed  is  that  of  the  10  parameters  measured,  the  value 
of  e(c)  ranks  fourth  in  significance. 


A  similar  analysis  regarding  risk  can  be  made  by 
computing  the  Significance  Ratios  for  the  Coefficient  of 
Variation.   In  this  case  Equation  1  can  be  rewritten  as: 

A  CV 


SR  = 


A  parameter 


(3) 


Table  12-12  contains  the  results  of  these  calculations. 

TABLE  12-12 
SIGNIFICANCE  RATIOS  -  CV 

Significance  Ratios 


Parameter 

Set 

Range 

Mean 

Normalized  Mean 

n 

.28 

-  1.36 

.67 

100 

P 

.31 

-  1.55 

.64 

96 

Q  at  5% 

.04 

-  .43 

.19 

28 

Q  at  951 

.01 

-  .18 

.06 

9 

Q (max) 

.06 

-  .15 

.10 

15 

Q(min) 

0 

0 

0 

C 
u 

.04 

-  .05 

.04 

6 

e(a) 

.06 

-  .11 

.09 

13 

e(b) 

0 

-  .06 

.03 

5 

e(c) 

.04 

-  .08 

.06 

9 
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From  Table  12-12  certain  observations  regarding  the 
significance  of  a  change  in  model  parameters  on  the  relative 
risk  of  the  venture  can  be  made.   First,  it  is  observed  that  a 
change  in  the  number  of  wells  drilled  and  in  the  probability 
of  success  both  have  about  the  same  effect  on  the  risk  of  the 
venture.   It  is  interesting  to  note  that  the  Normalized 
Significance  Ratios  for  E(P)  are  exactly  reversed  for  CV.   In 
this  regard,  Table  12-11  indicated  that  a  change  in  Q  at  951 
had  a  greater  significance  than  a  change  in  Q  at  5%.   Table 
12-12  reveals  the  opposite  is  true  for  the  relative  risk. 
Here  a  change  in  Q  at  5%    has  a  much  greater  significance  than 
a  change  in  Q  at  95%.   The  same  reversal  is  evident  in  the 
values  of  lowest  and  highest  profit  per  barrel,  e(a)  and  e(c). 

It  must  be  remembered,  of  course,  that  the  Signifi- 
cance Ratio  for  risk  (coefficient  of  Variation)  does  not  in- 
dicate whether  the  risk  is  increasing  or  decreasing  with  a 
parameter  change,  but  only  the  relative  amounts  that  it  is 
changing.   For  the  increase  or  decrease  correlation  we  must 
go  to  Tables  12-1  through  12-10  and  12-1A  through  12-9A. 

12.4   Convergence  of  the  Expected  Profit 

If  the  Monte  Carlo  simulation  model  is  executed  only 
once,  only  one  value  for  the  discounted  profit  is  generated  and 
there  is  no  distribution  of  values  that  can  be  analyzed.   Since 
simulation  is  being  used  because  the  model  is  too  complex  for 
other  types  of  analyses,  there  is  no  basis  for  relying  on  the 
single  result  as  a  typical  output  of  the  model  -  there  is  no 
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way  to  tell  whether  the  result  is  an  extreme  case  or  about  the 

average . 

As  the  simulation  process  is  repeated  frequencies  of 

occurrence  are  obtained  from  which  the  averages,  the  extremes, 

and  what  the  shape  of  the  distribution  would  be  if  the  model 

had  been  executed  infinitely  many  times  can  be  inferred.   In 

this  regard,  R.F.  Barton  [6]  poses  the  following  question: 

....  how  reliable  are  the  results  from  a  limited 
number  of  executions?   Two  executions  may  give  us  a 
great  deal  more  information  than  a  single  execution. 
One  hundred  iterations  may  give  us  hints  of  the  extre- 
mes and  shapes  of  the  ultimate  distributions.   A  thou- 
sand executions  may  perhaps  give  us  results  that  app- 
roximate smooth  curves.   But  a  second  thousand  execu- 
tions may  suggest  smooth  curves  of  different  locations 
and  shapes.   A  third  thousand  executions  may  generate 
results  that  look  even  different.   On  the  other  hand, 
the  separate  sets  of  one  thousand  executions  may  look 
so  similar  we  feel  no  need  for  further  executions. 

These  possibilities  raise  a  very  practical  question: 
How  long  to  run  a  simulation  model?   One  aspect  of  this  quest- 
ion was  analyzed  by  measuring  the  convergence  of  the  discounted 
Expected  Profit, E(P),  averaged  over  the  number  of  iterations 
run,  with  the  discounted  Expected  Profit  averaged  over  1000 
iterations . 

As  an  example,  if  the  value  E(P)  for  the  first  iteration 
is  $100.0  million  and  $50.0  million  for  the  second  iteration 
the  discounted  E(P)  as  determined  by  two  iterations  is  $75.0 
million.   We  can  then  compare  this  value  with  the  E(P)  averaged 
over  1000  iterations,  say  $60.0  million,  to  determine  the 
rapidity  of  convergence.   If  we  were  satisfied  with  an  Expected 
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Profit  that  varied  no  more  than  +  51  from  the  value  which  would 
be  obtained  with  1000  iterations  (90%  certainty)  then  by  plot- 
ting or  analyzing  the  data  we  could  determine  how  many  iter- 
ations were  necessary  to  reach  this  point.   It  may  turn  out,  in 
our  example,  that  a  value  of  $60.0  million  ±  $3.0  million  was 
achieved  after  only  500  iterations.   By  performing  this  analysis 
with  many  different  values  of  model  variables,  using  the  same 
random  numbers,  and  each  time  determining  the  point  beyond 
which  the  averaged  E(P)  varies  by  no  more  than  ±  51  of  the  1000 
iteration  value,  we  may  be  able  to  infer  the  number  of  itera- 
tions necessary  for  90%  confidence  limits  for  this  particular 
model . 

This  is,  in  fact,  what  has  been  done.   For  each  set  of 
model  variables  the  iteration  beyond  which  the  value  of  E(P) 
does  not  vary  by  more  than  ±  5%  from  the  value  of  discounted 
expected  profit  after  1000  iterations  was  measured.   These 
values  are  contained  on  Table  12-13.   Figures  12-9  through 
12-14  graphically  depict  this  convergence  process  for  various 
combinations  of  model  parameters. 

In  addition  to  this  analysis,  the  model  was  run  with 
the  base  data  for  a  total  of  10,000  continuous  iterations  to 
determine  if  any  periodicities  in  the  variables,  as  suggested 
by  Barton  [7] ,  occur  which  would  materially  effect  the  value 
of  the  discounted  Expected  Profit  after  10,000  iterations. 
Table  12-10  and  Figure  12-15  contain  the  results  of  this  test. 
It  can  be  seen  from  Figure  12-15  that  no  periodicities  occur 
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and  the  10,000  iteration  value  of  E(P)  is,  in  fact,  only 
slightly  different  from  the  1000  iteration  value. 

A  visual  inspection  of  the  data  contained  on  Table 
12-13  indicated  that  Expected  Profit,  E(P),  does  not  converge 
with  any  regularity,  either  within  80%  or  901  confidence  limits 
It  is  seen  that  even  within  a  parameter  set (i.e.  number  of 
wells  drilled)  the  number  of  iterations  vary  greatly.   The  only 
significant  fact  which  can  be  discerned  from  the  data  is  that 
values  of  E(P)  within  80%  confidence  limits  can  be  achieved 
anywhere  from  10  to  339  iterations  and  within  90%  confidence 
limits  anywhere  from  63  to  710  iterations.   If  the  highest 
and  lowest  iteration  value  is  eliminated  the  range  for  80% 
confidence  limits  is  18  to  242  iterations  and  the  range  for 
90%  confidence  limits  is  68  to  596  iterations. 
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Factor 

E(P) 

Case 

Being 

1000 

No.  of 

No. 

Varied 

Iterations 

Iterations 

TABLE  12-13 

CONVERGENCE  OF  E(P)  WITH  NUMBER  OF  ITERATIONS 
WITHIN  CONFIDENCE  LIMITS 


Within  80%  Conf.   Within  901  Conf 
Limits  (±  10%)     Limits  (±  5%) 

No.  of 
Iterations 


Base    -         68.02  117  207 

502 
457 
440 
268 
195 
107 
261 
238 
63 

710 

268 

395 

255 

83 

167 

99 

68 

96 

543 
596 
206 
207 
129 
129 
129 
129 

129 
207 
128 
207 
129 
207 
207 
473 


1 

n 

15.80 

239 

2 

n 

33.71 

242 

3 

n 

39.49 

190 

4 

n 

51.48 

186 

5 

n 

60.02 

147 

6 

n 

73.52 

65 

7 

n 

83.12 

111 

8 

n 

99.47 

136 

9 

n 

115.37 

10 

10 

P 

18.19 

339 

11 

p 

33.75 

213 

12 

p 

43.56 

204 

13 

P 

81.08 

108 

14 

P 

90.46 

41 

15 

P 

100.50 

56 

16 

P 

114.13 

42 

17 

p 

125.12 

26 

18 

P 

137.27 

40 

19 

Q(5%D 

46.15 

120 

20 

Q(5%) 

61.16 

117 

21 

Q(5%) 

63.24 

117 

22 

Q(5%) 

65.30 

117 

23 

Q(5%D 

69.34 

17 

24 

Q(5%) 

70.90 

17 

25 

Q(5%D 

72.89 

17 

26 

QC5I) 

80.42 

42 

27 

Q(95%) 

23.45 

105 

28 

Q(95%) 

52.42 

17 

29 

Q(95%) 

57.50 

17 

30 

Q(95%) 

62.64 

18 

31 

Q(95$) 

72.63 

117 

32 

Q(95%) 

77.14 

117 

33 

Q(95l) 

81.69 

118 

34 

Q(9  5%) 

97.46 

119 
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TABLE  12-13  Continued 


Within 

801  Conf. 

Within  90%  Conf. 

Limits     01) 

Limits  (±  51) 

Factor 

E(P) 

Case 

Being 

1000 

No.  of 

No.  of 

No. 

Varied 

Iterations 

Iterations 

Iterations 

35 

Q  max 

64.35 

43 

129 

36 

Q  max 

69.29 

18 

305 

37 

Q  min 

68.02 

117 

207 

38 

Q  min 

68.02 

117 

207 

39 

Q  min 

68.18 

22 

129 

40 

Cu 

69.30 

18 

129 

41 

Cu 

68.87 

18 

129 

42 

Cu 

67.17 

117 

207 

43 

Cu 

65.47 

118 

207 

44 

Cu 

62.07 

117 

265 

45 

Cu 

53.58 

119 

539 

46 

e  (a) 

60.99 

117 

207 

47 

e  (a) 

63.80 

117 

207 

48 

e  (a) 

65.21 

117 

129 

49 

e  (a) 

66.61 

117 

207 

50 

e  (a) 

69.43 

18 

207 

51 

e  (a) 

70.83 

19 

129 

52 

e  (a) 

72.24 

18 

129 

53 

e(b) 

60.74 

117 

129 

54 

e(b) 

62.97 

22 

129 

55 

e(b) 

65.50 

22 

129 

56 

e(b) 

70.54 

117 

207 

57 

e(b) 

73.07 

117 

207 

58 

e(b) 

75.31 

117 

126 

59 

e  (c) 

60.54 

117 

207 

60 

e  (c) 

61.65 

117 

207 

61 

e  (c) 

64.70 

117 

207 

62 

e  (c) 

71.35 

18 

129 

63 

e  (c) 

74.39 

22 

129 

64 

e  (c) 

77.71 

22 

129 

65 

e  (c) 

99.87 

22 

129 

Beta   Base 


67.98 


166 


239 
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12.5   Changes  in  Return  and  Risk  with  the  Number  of  Wells 
Drilled 

Since  the  risk  of  the  venture  has  been  measured  by  the 
coefficient  of  the  variation  it  may  be  of  interest  to  ascertain 
whether  there  is  a  mathematical  relationship  between  the  risk 
and  the  number  of  wells  drilled.   To  determine  whether  a  mathe- 
matical relationship  exists  we  start  with  Equation  7,  Chapter  V 


CV  =  i  Ji-£ 

Now,  if  p  is  held  constant  it  can  be  seen  that  the  Coefficient 
of  Variation  should  be  proportional  to  the  number  of  wells 
drilled  as  follows: 

CV  ~  -i  (4) 

/n 

or, 

CV  x  /n  =  Constant  (5) 

Table  12-14  contains  a  test  of  this  proportionality  hypothesis 
using  the  data  contained  on  Table  12-1 
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TABLE  12-14 
PROPORTIONALITY  TEST  -  CV  to  n 

n  /n      Coeff.  of  Variation     Constant 

1.59  3.54 

1.11  3.52 

1.07  3.70 

0.92  3.56 

0.84  3.56 

0.79  3.53 

0.75  3.52 

0.71  3.55 

0.64  3.50 

0.62  3.67 

It  can  be  seen  from  Table  12-14  that  the  Coefficient  of  Vari- 
ation is  indeed  proportional  to  1//H.   Thus,  if  the  number  of 
wells  is  increased  by  4  times  the  risk  will  be  reduced  by  a 
factor  of  2. 

It  is  also  interesting  to  note  that  the  Best  Estimate 
of  Discounted  Profit  per  well,  or  E(P)/n,  remains  fairly  con- 
stant while  the  risk,  as  measured  by  the  Coefficient  of  Vari- 
ation, decreases  with  l//n.   This  can  be  seen  on  Table  12-15. 
This  would  tend  to  indicate  that  while  the  best  estimate  of 
return  per  well  remains  fairly  constant  with  the  number  of 
wells  drilled  the  risk  associated  with  that  return  deminishes 
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as  we  increase  the  number  drilled. 

TABLE  12-15 
BEST  ESTIMATE  OF  PROFIT  PER  WELL  AND  RISK 


Number  of    Best  Estimate        Coeff.  of   Constant 
Wells  Drilled   of  Discounted   E (P)   Variation     1 
Case       (n)        Profit-  E(P)    n        CV       /^ 


15.80  3.16  1.59  3.54 

33.71  3.37  1.11  3.52 

39.49  3.29  1.07  3.70 

51.48  3.43  0.92  3.56 

60.02  3.33  0.84  3.56 

68.02  3.40  0.79  3.53 

73.52  3.34  0.75  3.52 

83.12  3.32  0.71  3.55 

99.47  3.32  0.64  3.50 

115.37  3.29  0.62  3.67 
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CHAPTER  XIII 
CONCLUSIONS 

13.1   General 

The  following  general  conclusions  are  drawn  from  the 
results  of  this  study: 

1.  The  model  developed  herein  provides  a  suitable 
vehicle  for  measuring  the  sensitivity  of  the  parameters  affect- 
ing the  profitability  of  an  exploratory  oil  drilling  venture. 

2.  The  parameters  in  the  profitability  model  should 
include  not  only  single  best  estimate  values  but,  where  app- 
ropriate, also  a  probability  distribution  of  values. 

3.  The  probability  distribution  best  representing  the 
number  of  successes  is  the  Binomial  distribution,  especially 
when  the  probability  of  success  and  number  of  wells  drilled  are 
in  the  order  of  magnitude  normally  encountered  in  an  exploratory 
drilling  program.   For  probabilities  of  success  approaching  501 
and  large  numbers  of  wells  drilled,  the  Normal  or  Gaussian  dis- 
tribution can  closely  approximate  the  Binomial  distribution. 

4.  The  reserve  size  discovered  by  a  successful  well  is 
best  represented  by  a  Lognormal  distribution.   The  values  can 
be  calculated  manually  using  points  scaled  from  a  straight  line 
drawn  on  log-probability  paper  between  values  of  reserve  size 
which  can  be  expected  with  5  and  951  probability.   The  lognormal 
distribution  can  be  simulated  in  computer  applications  with 
excellent  correlation  using  the  Lognormally  Distributed  Variable 
Generator  contained  in  this  study. 
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5.  The  profit  per  barrel  received  over  a  period  of 
years  from  oil  produced  must  be  discounted  to  reflect  the  time 
value  of  money. 

6.  The  Triangular  distribution  function  adequately 
serves  to  represent  the  discounted  profit  per  barrel.   The  Beta 
distribution  function  provides  only  slightly  different  values 
of  discounted  profit  per  barrel  when  tested  using  the  same 
input  data  used  in  the  Triangular  distribution.   Both  the 
Triangular  and  Beta  distributions  can  be  simulated  in  computer 
applications  using  the  Generators  contained  herein. 

7.  The  complexity  of  the  distributions  in  the  model 
precludes  manual  computation  of  profitability  and  sensitivity 
studies  thereby  necessitating  the  use  of  the  Monte  Carlo 
technique  and  a  digital  computer. 

13.2   Sensitivity  Studies 

In  the  sensitivity  studies,  where  individual  parameters 
in  the  model  were  varied,  with  the  remaining  parameters  being 
held  constant,  the  following  conclusions  were  drawn: 

1.  The  one  single  observed  value  which  fluctuated  the 
greatest  (on  a  percentage  basis)  was  the  minimum  value  of 
discounted  profit  within  90%  Confidence  Limits. 

2.  Except  in  isolated  cases,  the  chance  if  a  profit 
greater  than  zero  in  all  tests  remained  fairly  constant  at  95 
to  96%.   Conversely,  there  remained  a  fairly  constant  chance 
of  4  to  5%  of  incurring  a  loss. 
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3.  A  subsitiution  of  the  Beta  distribution  for  dis- 
counted profit  per  barrel  rather  than  the  Triangular  distri- 
bution makes  very  little  effect  on  the  best  estimate  of  dis- 
counted profit  for  the  venture  or  the  risk  as  measured  by  the 
coefficient  of  variation. 

4.  There  appears  to  be  no  periodicities  in  the  para- 
meters of  the  model  and  the  best  estimate  of  discounted  profit 
calculated  after  10,000  executions  varies  only  slightly  from 
the  value  calculated  after  1,000  iterations. 

5.  It  is  possible  to  calculate  values  of  the  best 
estimate  of  discounted  profit  for  the  venture  within  801 
confidence  limits  (±  101  of  the  1,000  iteration  value  of  E(P)) 
by  executing  the  model  between  18  and  242  times;  and  within 
90%  confidence  limits  by  executing  the  model  between  68  and 
596  iterations.   Further,  there  appears  to  be  no  regular 
convergence  of  E(P)  with  the  number  of  iterations  performed. 

6.  A  Significance  Ratio,  relating  a  change  in  observed 
value  calculated  by  the  model  to  a  change  in  parameter  value, 
provides  an  excellent  tool  for  comparing  the  effects  of  changes 
of  different  model  parameters,  especially  when  the  Significance 
Ratio  is  normalized  on  a  1  to  100  scale. 

7.  Using  the  Significance  Ratios  to  relate  parameter 
changes  to  discounted  venture  profitability,  E(P)  it  is  con- 
cluded that  with  the  parameter  values  used: 

a.   Changes  in  probability  of  success  and  number 
of  wells  drilled  have  almost  equal  significance  on  changes  in 
venture  profitability. 
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b.  Changes  in  the  lower  range  of  reserve  size  have 
only  about  30$  of  the  effect  of  changes  in  p  or  n  on  E(P). 
However,  changes  in  the  upper  range  have  a  70%  significance. 
The  values  selected  at  lower  and  upper  boundaries  for  reserve 
sizes  in  the  model  have  little  or  no  effect  on  the  value  or  E(P) 

c.  Changes  in  the  highest  value  of  discounted 
profit  per  barrel  have  more  than  twice  the  significance  as 
changes  in  the  lowest  values  of  discounted  profit  per  barrel; 
but  even  then,  have  only  461  of  the  significance  as  changes 
in  n  or  p  on  E (P) . 

d.  Changes  in  dry  hole  cost  have  relatively 
little  effect  on  E(P) . 

8.   Using  the  Significance  Ratios  to  relate  parameter 
changes  to  the  relative  risk  of  the  venture  as  measured  by  the 
coefficient  of  variation,  CV,  it  is  concluded  that  with  the 
parameter  values  used: 

a.  Changes  in  the  number  of  wells  drilled,  n,  and 
the  probability  of  success,  p,  have  about  the  same  effect  on 
the  risk  of  the  venture. 

b.  Changes  in  the  lower  range  of  reserve  size  are 
next  in  significance.   Changes  in  this  parameter  are  more  than 
3  times  as  significant  as  changes  in  the  upper  ranges  of 
reserve  sizes. 

c.  The  value  selected  as  the  upper  boundary  for 
reserve  size  will  effect  the  risk  in  the  venture  but  the  lower 
boundary  limit  will  not. 
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d.  Changes  in  dry  hole  costs  do  little  to  change 
the  risk  of  the  venture. 

e.  Changes  in  the  lowest  value  of  discounted 
profit  per  barrel  is  about  1.5  times  as  significant  as  changes 
in  the  highest  value  of  discounted  profit  per  barrel  but  even 
then,  has  only  a  significance  of  131  as  compared  with  the  risk 
associated  with  changes  in  the  number  of  wells  drilled. 

9.   While  the  return  per  well  remains  fairly  constant 
as  the  number  of  wells  increase,  the  relative  risk  of  the 
venture  decreases  as  the  square  root  of  the  number  of  wells 
drilled. 

The  seemingly  large  number  of  findings  drawn  from  the 
results  of  this  study  are  indicative  of  the  value  of  this  type 
of  analysis.   It  is  recognized,  of  course,  that  not  all  the 
conclusions  listed  herein  will  be  considered  of  equal  import- 
ance by  the  reader  -  it  would  be  rare  if  they  did.   The  fact 
that  these  findings  could  be  drawn  at  all,  however,  leads  to 
the  overall  conclusion  that  sensitivity  analysis  of  this  type 
provides  a  valuable  tool  for  the  decision-maker. 


CHAPTER  XIV 
RECOMMENDATIONS 

It  is  in  the  light  of  the  problems  facing  most  decision- 
makers, namely  the  "What  if...?"  questions  cited  in  Chapter  XI 
and  the  contention  so  lucidly  expressed  by  William  T.  Morris 
[35]  that  at  some  point  we  must  stop  gathering  data  and  re- 
fining our  information  and  start  recommending  action,   that  the 
profitability  model  is  ideally  suited. 

It  is  therefore  recommended  that  a  model  of  this  type 
be  used  not  to  obtain  absolute  data  on  whether  a  profit  will 
be  made  for  the  venture  (or  some  other  observation  of  this 
type)  since  other  methods  of  analysis  may  be  better  suited  for 
this  purpose.   Instead,  it  is  recommended  that  the  model  be 
used  to  assess  the  relative  sensitivity  of  the  results  to  para- 
meter changes  in  the  model.   With  this  knowledge  of  parameter 
sensitivity  we  can  now  devote  time  and  effort  to  refining  only 
those  parameters  which  have  a  significant  affect  on  the  outcome. 

For  example,  if  we  were  concerned  with  the  best  estimate 
of  discounted  profit  for  the  venture  we  would  probably  devote 
our  efforts  to  refining  our  estimates  of  the  probability  of 
success,  the  upper  range  of  reserve  size  and  the  highest  value 
of  discounted  profit  per  barrel  (assuming,  of  course,  that  the 
number  of  wells  drilled  is  fixed).   We  would  not,  in  this 
example,  concern  ourselves  with  such  items  as  the  dry  hole 
costs,  the  lower  range  of  reserve  sizes  and  the  like. 
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In  this  way,  using  the  sensitivity  analysis  approach, 
the  decision-maker  can  achieve  the  maximum  return  for  the 
effort  he  can  afford  to  expend. 


APPENDIX  A 
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THE  LOGNORMALLY  DISTRIBUTED  VARIABLE  GENERATOR 

The  FORTRAN  IV  Computer  program  for  lognormally 
distributed  deviates  requires  two  different,  uniformally 
distributed  random  fractions  to  generate  one  normal  variable 
With  two  uniformally  distributed  fractions,  RA  and  RB,  one 
generates  a  variable,  V,  which  is  normally  distributed  about 
a  mean  of  0,  with  a  standard  deviation  of  1,  thus, 

V  =  (-2  In  RA)  '5  x  COS  2tt  (RB) 

To  convert  V  to  a  variable  QNORM,  which  is  normally  distrib- 
uted about  a  mean  of  TMU,  with  a  standard  deviation  SD,  we 
simply  use  the  relationship, 

QNORM  =  (V)  SD  +  TMU 

By  multiplying  the  natural  logarithm  base,  e,  to  the  normal 
deviate  power,  a  lognormally  distributed  deviate,  QLOGN,  is 
generated. 


QLOGN  =  e 


QNORM 


Once  the  Lognormally  distributed  deviate  is  generated 
it  is  tested  against  the  lower  and  upper  limils  established  as 
inputs  to  the  program.   If  the  deviate  is  less  than  the  minimum 
value  of  QMIN,  it  is  set  to  read  QMIN.   If  it  is  greater  than 
the  maximum  value  of  QMAX,  it  is  set  to  read  QMAX. 

It  will  be  noted  that  the  program  requires  a  value  of 
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mean  and  standard  deviation.   These  values  can  be  obtained 
from  the  median  value  of  Q  at  W(  ^Q)  ■  50%,  and  the  value  of 
Q  at  W(<Q)  =  84%  (or,  ym  +  a   )    as  follows: 

From  Equation  3, 

ym  (median)  =  yQ  (mean) 
From  Equation  4B, 

ln  xm  =  y0 

Qmedian  =  xm 
and,         TMU  =  y0 


then, 


TMU  =  loge(Qmedian) 
or 

TMU  =  ALOG (QMEDIAN) 

^o   ^m 


^m  +  °y   =  ^o  +  °y 


ym  +  ay  =  ln  Q  at  W(<Q)  =  84% 
or, 

ym  +  ay   =  ln  Q(84) 
and, 

ln  xm  =  ^m 
Therefore, 

^m  +  V  "  Xm  =  ln  Q(84)  "  lnCQm) 
or, 

SD  =  ALOG(Q84)  -  ALOG (QMEDIAN) 
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The  FORTRAN  IV  Computer  program  incorporating  these  concepts 
is  listed  as  follows: 

PROGRAM  QLOGN  (INPUT , OUTPUT) 

READ  100,  QMEDIAN,  Q84,  QMAX,  QMIN 

SD  =  ALOG(Q84)  -  ALOG (QMEDIAN) 

TMU  =  ALOG (QMEDIAN) 

RA  =  RANF  (0) 

RB  =  RANF  (0) 

V  =  (-2.0*ALOG(RA)**.5*  COS  (2 . 0*3 . 1416*RB) 

QNORM  =  V*SD  +  TMU 

QLOGN  =  EXP (QNORM) 

IF  (QLOGN  -  QMIN)  6,  7,  8 

6  QLOGN  =  QMIN 
GO  TO  7 

8  IF  (QLOGN  -  QMAX)  7,  7,  9 

9  QLOGN  =  QMAX 

7  CONTINUE 

PRINT  120,  QLOGN 
100   FORMAT  (10X,4F10.3) 
120   FORMAT  (20X,F10.3) 

END 
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THE  TRIANGULAR  DISTRIBUTION  VARIABLE  GENERATOR 

The  FORTRAN  IV  Computer  program  for  deviates  with  a 
Triangular  distribution  requires  one  uniformally  distributed 
random  function,  RANF(O).   This  uniformally  distributed  random 
function  will  assume  either  the  value  of  W(x  <  x^) ,  -  WXLTXI  - 
or  W(x  >  x-)  -  WXGTXI  -  depending  on  whether  the  calculated 
value  of  x-(XI)  is  equal  to  or  less  than  the  inputed  value  of  b 

The  FORTRAN  IV  program  incorporating  these  concepts  is 
listed  as  follows: 

PROGRAM  TRIA  (INPUT, OUTPUT) 

THIS  PROGRAM  CALCULATES  THE  BEST  ESTIMATE  OF  PROFIT/BBL 
ASSUMING  A  *  TRIANGULAR  *  DISTRIBUTION  OF  PROBABILITY. 

A  =  MINIMUM  PROFIT/BBL  EXPECTED 

B  =  MOST  LIKELY  VALUE  OF  PROFIT/BBL  EXPECTED 

C  =  MAXIMUM  PROFIT/BBL  EXPECTED 

READ  112,  A,  B,  C, 

CALL  RANF(O) 

WXLTXI  =  RANF(O) 

WXGTXI  =  1.0  -  WXLTXI 

XI  =  A  +  SQRT(WXLTXI*(B-A)*(C-A)) 

IF  (XI.LE.B)  55,50 
50   XI  =  C  -  SQRT(WXGTXI*(C-B)*C-A)) 
55   PROFIT  =  XI/100. 

PRINT  300,  PROFIT 
112   FORMAT  (10X,3F10.2) 
300   FORMAT  (15X,F4.2) 

END 
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THE  BETA  DISTRIBUTION  VARIABLE  GENERATOR 

The  FORTRAN  IV  Computer  program  for  deviates  with  a 
Beta  distribution  is  listed  as  follows: 


PROGRAM  BETA  (INPUT , OUTPUT) 

THIS  PROGRAM  CALCULATES  THE  BEST  ESTIMATE  OF  PROFIT/BBL 
ASSUMING  A  *  BETA  *  DISTRIBUTION  OF  PROBABILITY. 

A  =  MINIMUM  PROFIT/BBL  EXPECTED 

B  =  MOST  LIKELY  VALUE  OF  PROFIT/BBL  EXPECTED 

C  =  MAXIMUM  PROFIT/BBL  EXPECTED 

READ  100,  A,  B,  C 

XMU  =  (4.*B+A+C)/6.0 

XVAR  =  (C-A)**2/36.0 

BMEAN  =  (XMU-A)/(C-A) 

BVAR  =  XVAR/(C-A)**2 

XK1  =  BMEAN* (BMEAN* (1.0 -BMEAN) /BVAR- 1.0) 

XK2  =  XK1*((1.0-BMEAN) /BMEAN) 

CALL  GAMMARN  (XK1,GAM1) 

CALL  GAMMARN  (XK2,GAM2) 

BETA  =  (GAM1/(GAM1  +  GAM2))*(C-A)  +A 

PROFIT  =  BETA/100  . 

PRINT  300,  PROFIT 
100   FORMAT  (10X,3F10.2) 
300   FORMAT  (15X,F4.2) 

END 

SUBROUTINE  GAMMARN (TK, GAM) 

GAMMA  =1.0 

Kl  =  TK 

TK1  =  Kl 

Yl  =  RANF(O) 

IF  (Yl- (TK-TK1))  10,  10,  20 
10   Kl  =  Kl  +  1 
20   DO  30  I  =  1,K1 

Y2  =  RANF(0) 
30   GAMMA  =  GAMMA* Y 2 

GAM  =  -ALOG (GAMMA) 

RETURN 

END 
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COMPUTER  PRINTOUT  FOR  BASE  DATA 


PROFIT  DISTRIBUTION  OF  AN  EXPLORATORY  DRILLING  PROGRAM 


NO.  OF  WELLS    PROH  OF   DRY  HOLE 

DRILLED      SUCCESS     COST     QMEDIAN    Q84      QMAX      QMIN 
20         .150     .150      20.00     46.00    140.00      1.00 


THE  TRIANGULAR  DISTRIBUTION  FOR  DETERMINATION  OF  THE 
PROFIT  PER  BARREL  HAS  THE  FOLLOWING  EXTREMES  IN  CENTS  - 
A  =   50.00        B  =   88.00        C  =  115.00 


ITERATION 

NUMBER  OF 

FIELD 

PROFIT  PER 

ITERATION 

BEST  ESTIMATE 

SUCCESSES 

SIZES 

BARREL 

PROFIT 

OF  PROFIT 

1 

2 

52.90 

.83 

41.17 

41.17 

2 

3 

64.03 

1.01 

61.90 

51.53 

3 

0 

0.00 

0.00 

-3.00 

33.36 

4 

4 

114.52 

1.02 

114.13 

53.55 

5 

5 

189.57 

.79 

148.40 

72.52 

6 

0 

0.00 

0.00 

-3.00 

59.93 

7 

3 

34.12 

.87 

27.05 

55.24 

8 

5 

254.88 

1.01 

254.59 

80.15 

9 

3 

198.18 

•  81 

157.93 

88.80 

10 

2 

32.61 

1.01 

30.36 

82.95 

11 

0 

0.00 

0.00 

-3.00 

75.14 

12 

2 

38.72 

.93 

33.17 

71.64 

13 

5 

179.08 

.93 

164.57 

78.79 

14 

3 

95.99 

.57 

51.88 

76.87 

15 

3 

166.86 

.61 

98.92 

78.34 

16 

3 

129.04 

.65 

81.67 

78.55 

17 

2 

28.20 

•  86 

21.46 

75.19 

18 

3 

54.65 

•  80 

41.00 

73.29 

19 

3 

73.73 

.95 

67.13 

72.97 

20 

1 

85.62 

.97 

80.58 

73.35 

21 

5 

127.41 

.81 

100.96 

74.66 

22 

1 

3.84 

.92 

.69 

71.30 

23 

2 

40.45 

.79 

29.10 

69.46 

24 

4 

59.44 

.87 

49.05 

68.61 

25 

2 

53.91 

•  81 

41.17 

67.52 

26 

3 

49.27 

•  66 

30.22 

66.08 

27 

2 

26.85 

•  70 

16.03 

64.23 

28 

5 

96.30 

.73 

67.94 

64.36 

29 

2 

27.60 

.74 

17.80 

62.75 

30 

4 

187.40 

.99 

183.55 

66.78 

31 

2 

137.95 

.64 

85.37 

67.38 

32 

4 

83.15 

.82 

65.70 

67.33 

33 

2 

107.08 

.82 

84.63 

67.85 

34 

0 

0.00 

0.00 

-3.00 

65.77 

35 

4 

84.91 

.85 

69.90 

65.89 

36 

2 

38.40 

.94 

33.52 

64.99 

37 

2 

25.44 

1.02 

23.35 

63.86 

38 

4 

125.78 

.73 

89.62 

64.54 

39 

1 

6.76 

.78 

2.41 

62.95 

40 

5 

110.73 

.89 

95.83 

63.77 

41 

0 

0.00 

0.00 

-3.00 

62.14 

42 

4 

116.96 

1.09 

124.90 

63.63 

43 
44 
45 
46 
47 
48 

49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 

70 
71 
72 

73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 
91 
92 
93 
94 
95 
96 
97 
98 
99 
100 
101 
102 

103 
104 


6 

229.75 

.91 

206.55 

66.96 

1 

32.97 

.97 

29.02 

66.10 

2 

25.19 

1.07 

24.20 

65.17 

0 

0.00 

0.00 

-3.00 

63.68 

2 

181.65 

1.02 

183.29 

66.23 

3 

84.76 

.87 

70.92 

66.33 

5 

138.59 

.59 

79.57 

66.60 

8 

160.84 

1.04 

165.66 

68.58 

4 

91.12 

1.04 

92.52 

69.05 

1 

50.01 

.71 

32.64 

68.35 

2 

30.00 

.67 

17.45 

67.39 

1 

43.41 

.55 

20.90 

66.53 

5 

129.75 

.94 

119.35 

67.49 

1 

19.93 

.79 

12.84 

66.51 

4 

63.36 

.92 

55.84 

66.32 

3 

72.45 

.87 

60.17 

66.22 

3 

110.63 

.63 

66.80 

66.23 

3 

100.93 

.66 

64.42 

66.20 

5 

194.06 

.94 

181.04 

68.08 

3 

170.74 

.70 

117.80 

68.88 

0 

0.00 

0.00 

-3.00 

67.74 

3 

74.02 

•  86 

60.82 

67.63 

2 

39.45 

.87 

31.43 

67.08 

4 

89.92 

.96 

83.53 

67.32 

3 

112.67 

.84 

91.91 

67.69 

6 

172.22 

•  88 

150.10 

68.90 

3 

51.86 

.99 

48.86 

68.61 

2 

47.19 

.84 

36.78 

68.16 

3 

63.76 

.91 

55.78 

67.98 

2 

76.24 

.58 

41.86 

67.62 

2 

23.86 

•  80 

16.45 

66.92 

3 

110.39 

•  80 

85.40 

67.17 

1 

20.76 

.84 

14.53 

66.47 

3 

95.34 

.93 

85.91 

66.72 

2 

38.74 

.87 

30.88 

66.26 

4 

23.97 

.76 

15.85 

65.61 

3 

32.79 

.79 

23.31 

65.08 

3 

50.87 

.89 

42.79 

64.80 

2 

91.19 

.73 

64.26 

64.79 

5 

96.51 

•  63 

58.39 

64.71 

5 

119.23 

•  91 

106.21 

65.21 

0 

0.00 

0.00 

-3.00 

64.40 

5 

184.12 

.73 

132.41 

65.20 

2 

31.18 

.78 

21.46 

64.69 

3 

51.86 

•  86 

41.95 

64.43 

3 

80.81 

.58 

44.18 

64.20 

3 

109.49 

.95 

101.59 

64.62 

3 

55.49 

.57 

28.98 

64.23 

2 

19.07 

1.12 

18.70 

63.72 

2 

29.24 

.94 

24.86 

63.30 

3 

88.87 

.70 

59.78 

63.26 

1 

23.84 

•  88 

18.05 

62.78 

3 

83.49 

.67 

53.13 

62.68 

1 

88.52 

1.08 

92.78 

63.00 

4 

156.23 

•  84 

129.02 

63.68 

0 

0.00 

0.00 

-3.00 

63.00 

3 

50.68 

1*05 

50.79 

62.87 

2 

42.56 

.93 

37.08 

62.61 

0 

0.00 

0.00 

-3.00 

61.96 

3 

55.47 

.87 

45.67 

61.81 

4 

88.88 

.91 

78.40 

61.97 

1 

34.29 

.89 

27.59 

61.64 

105 
106 
107 
108 
109 
110 
111 
112 
113 
114 
115 
116 

117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 
130 
131 
132 
133 
134 
135 
136 
137 

133 
139 

140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
151 
152 

153 
154 
155 
156 
157 
158 
159 
160 
161 

162 
163 

164 

165 
166 


5 

124.19 

•  68 

82.68 

61.84 

4 

82.08 

1.13 

90.48 

62.11 

2 

92.60 

.65 

57.52 

62.06 

4 

141.82 

.75 

103.75 

62.45 

5 

101.35 

.87 

85.45 

62.66 

3 

45.26 

.77 

32.31 

62.38 

4 

76.50 

.77 

56.34 

62.33 

2 

27.12 

.99 

24.04 

61.99 

1 

18.92 

.96 

15.36 

61.58 

7 

83.56 

.83 

67.61 

61.63 

3 

31.66 

.84 

24.00 

61.30 

2 

85.46 

.64 

52.03 

61.22 

6 

166.72 

.84 

137.93 

61.88 

5 

199.84 

.68 

133.63 

62.49 

5 

118.33 

.94 

108.42 

62.87 

3 

150.97 

.85 

125.42 

63.39 

0 

0.00 

0.00 

-3.00 

62.84 

4 

256.55 

.60 

151.98 

63.57 

8 

199.11 

.93 

184.32 

64.56 

2 

111.75 

.84 

91.60 

64.77 

2 

25.56 

.80 

17.70 

64.40 

0 

0.00 

0.00 

-3.00 

63.86 

2 

57.89 

.73 

39.42 

63.67 

6 

83.49 

1.03 

83.96 

63.83 

5 

245.33 

.95 

229.78 

65.11 

3 

194.04 

.95 

181.68 

66.01 

5 

108.37 

1.02 

108.45 

66.34 

3 

72.57 

.54 

36.64 

66.11 

2 
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THE  BEST  ESTIMATE  OF  PROFIT  BY  THE  MONTE  CARLO  SIMULATION  METHOD 
AT  1000  ITERATIONS  IS    68.02  MILLION  DOLLARS 
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THE  VARIANCE  IS      2879.367 

THE  STANDARD  DEVIATION  IS        53.660 

THE  COEFFICIENT  OF  VARIATION  IS   .7889 


PROFIT  DISTRIBUTION  OF  AN  EXPLORATORY  DRILLING  PROGRAM 


NO.  OF  WELLS 

PROB  OF 

DRY  HOLE 

DRILLED 

SUCCESS 

COST 

QMEDIAN 

Q84 

QMAX 

QMIN 

20 

.150 

.150 

20.00 

46.00 

140.00 

1.00 

THE  TRIANGULAR  DISTRIBUTION  FOR  DETERMINATION  OF  THE 
PROFIT  PER  BARREL  HAS  THE  FOLLOWING  EXTREMES  IN  CENTS  - 
A  =   50.00        B  =   88.00        C  =  115.00 
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